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Abstract 

An increasing number of studies demonstrate phenotypic and genetic changes in 
natural populations that are subject to climate change, and there is hope that 
some of these changes will contribute to avoiding species extinctions ('evolution- 
ary rescue'). Here, we review theoretical models of rapid evolution in quantitative 
traits that can shed light on the potential for adaptation to a changing climate. 
Our focus is on quantitative-genetic models with selection for a moving pheno- 
typic optimum. We point out that there is no one-to-one relationship between 
the rate of adaptation and population survival, because the former depends on 
relative fitness and the latter on absolute fitness. Nevertheless, previous estimates 
that sustainable rates of genetically based change usually do not exceed 0.1 
haldanes (i.e., phenotypic standard deviations per generation) are probably cor- 
rect. Survival can be greatly facilitated by phenotypic plasticity, and heritable var- 
iation in plasticity can further speed up genetic evolution. Multivariate selection 
and genetic correlations are frequently assumed to constrain adaptation, but this 
is not necessarily the case and depends on the geometric relationship between the 
fitness landscape and the structure of genetic variation. Similar conclusions hold 
for adaptation to shifting spatial gradients. Recent models of adaptation in multi- 
species communities indicate that the potential for rapid evolution is strongly 
influenced by interspecific competition. 



Introduction 

Over the past two decades, it has become clear that evolu- 
tionary change can be fast enough to be observed in pres- 
ent-day populations (Hendry and Kinnison 1999; Kinnison 
and Hendry 2001; Hendry et al. 2008; Gingerich 2009) and 
that it can directly affect the dynamics of populations and 
communities (Hairston et al. 2005; Saccheri and Hanski 
2006; Kinnison and Hairston 2007; Pelletier et al. 2009). 
Much recent interest has focused on the possibility that 
so-called rapid or contemporary evolution leads to 'evolu- 
tionary rescue', whereby threatened populations avoid 
extinction by adapting to an altered environment (Barrett 
and Hendry 2012; Gonzalez et al. 2013). This issue is par- 
ticularly pressing in the context of global climate change, 
which subjects large numbers of populations to shifts in 
temperature, aridity, seasonal patterns, etc. While pheno- 
typic responses to climate change have been documented 



(Bradshaw and Holzapfel 2006; Parmesan 2006; Hoffmann 
and Sgro 2011 and this issue), the potential for evolution- 
ary rescue is still unclear (Bell 2013). At the same time, it is 
often difficult to distinguish changes based on genetic evo- 
lution from those due to phenotypic plasticity (Merila 
2012; Merila and Hendry 2014). 

At the basis of many questions in the context of adap- 
tation to environmental change are rates of phenotypic 
evolution (Hendry and Kinnison 1999; Kinnison and 
Hendry 2001; Gingerich 2009). These rates are often 
measured in haldanes. One haldane is equivalent to a 
change in one phenotypic standard deviation per genera- 
tion (for other measures, see discussion in Hendry and 
Kinnison 1999, and for alternative standardizations and 
issues of scale, Hereford et al. 2004; Hansen and Houle 
2008). Several recent meta-analyses of contemporary evo- 
lution yield the following picture: evolutionary rates 
above 0.1 haldanes are not uncommon (Hendry and 
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Kinnison 1999; Gingerich 2009), even though the major- 
ity of rates are lower (Kinnison and Hendry 2001). Rates 
are higher in populations that are strongly influenced by 
human activities (Hendry et al. 2008; Darimont et al. 
2009). Rates measured over few generations are higher 
than those measured over many generations (Gingerich 
1983; Kinnison and Hendry 2001; Hendry et al. 2008; 
Gingerich 2009; Westley 2011). Studies that controlled 
for environmental effects (e.g., using common garden 
experiments) find lower rates than those that do not 
(Hendry et al. 2008), suggesting a role for phenotypic 
plasticity (Pigliucci and Murren 2003; Hendry et al. 2008; 
Westley 2011). Over paleontological timescales, the best- 
fitting model of phenotypic evolution is one of stasis 
interrupted by bursts of change (Estes and Arnold 2007; 
Uyeda et al. 2011). 

The aim of this study is to review quantitative-genetic 
models that shed light on the potential for rapid adapta- 
tion. Our focus will be on the evolution of quantitative 
traits, that is, traits with continuous variation that are 
determined by a large number of loci with appreciable 
standing genetic variation. While we will frequently men- 
tion the link between adaptation and population survival, 
we do not aim for a comprehensive review of evolutionary 
rescue theory (see Gonzalez et al. 2013 and 14 other articles 
in a recent theme issue of the Philosophical Transactions of 
the Royal Society B, vol. 368:1610). In particular, we will 
not treat evolutionary rescue via the fixation of single large 
mutations (Gomulkiewicz and Holt 1995; Holt and Gom- 
ulkiewicz 1997; Orr and Unckless 2008; Uecker and Her- 
misson 2011; Kirkpatrick and Peischl 2013; Martin et al. 
2013). 

The structure of the article is as follows. We first give a 
detailed description of the basic models of adaptation of 
single and multiple quantitative traits under various sce- 
narios of environmental change, including a discussion of 
'maximal sustainable rates of evolution' (Burger and Lynch 
1995). Subsequently, we discuss four avenues into which 
the basic models have been extended by recent work (i) the 
role of phenotypic plasticity and its interactions with 
genetic evolution, (ii) determinants of adaptive potential 
and evolvability, (iii) adaptation to shifting spatial gradi- 
ents and (iv) evolution and adaptation in a community 
context. 

Basic models 

Modeling approaches 

Environmental change 

Most theoretical approaches to adaptation in a changing 
environment are based on models of stabilizing selection 
with a moving optimum. That is, at any given time, selec- 



tion favors a specific trait value (or combination of trait 
values), but this favored phenotype changes over time. The 
most important scenarios are the following: 

• A single, sudden change in the optimum: this is a classic 
scenario studied in population genetics and also in 
recent models about the genetic basis of adaptation 
(Orr 2005) and evolutionary rescue (Orr and Unckless 
2008). It is well suited to study adaptation in invasive 
species, as well as in species suffering a sudden degrada- 
tion of their environment. 

• Gradual (typically linear) movement of the optimum: 
this scenario seems best suited to investigate the effects 
of continued climate change (Fig. 1). 

• Random fluctuations of the optimum, either around a 
constant value or around a linear trend: these fluctua- 
tions may or may not show autocorrelation. Such mod- 




Figure 1 Illustration of trait evolution in the one-dimensional moving- 
optimum model. (A) Solid and dotted gray curves represent the fitness 
landscape at different points in time (eq. 2), whose width is determined 
by co. 8, is the optimal phenotype, which moves at constant speed 
(B t = kt). The black curves represent the distribution of breeding values 
in the population (mean z t , variance op. The mean phenotype evolves 
according to eqn (1). At the dynamic equilibrium, it follows the opti- 
mum with a constant lag S* t . (B) illustrates the relation between rate of 
evolution and extinction risk. The gray curves show the log mean fitness 
as a function of the mean phenotype z t for two different fitness func- 
tions with widths w s and w w , respectively. The rate of evolution, given 
by the horizontal arrows is determined by the fitness gradient fS t , indi- 
cated by the black lines. The vertical position of the population gives its 
mean (log) fitness. In the Figure, the optimum is assumed to move at 
rate k = 0.03B, and the population placed at the narrow fitness curve 
follows at this pace while maintaining a positive growth rate (w > 1). 
With the wide fitness function, however, the same rate of evolution 
requires a larger distance from the optimum, such that the growth rate 
is negative and the population goes extinct. 
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els are useful to study the effects of environmental sto- 
chasticity that overlay all climate-driven trends. 

Genetic adaptation 

The majority of models reviewed here are based on quanti- 
tative genetics theory. Evolving traits are assumed to have a 
polygenic basis and follow a normal distribution with phe- 
notypic variance a 2 . In the simplest case (additive genetics, 
no phenotypic plasticity), a 2 can be decomposed into 
a 2 = a 2 + a 2 , where a 2 is the additive genetic variance, a 2 
is the environmental variance (variation due to develop- 
mental instability and micro-environmental fluctuations), 
and h 2 = o 2 g /a 2 is the (narrow-sense) heritability. If pheno- 
types are measured in units of the environmental variance, 
a 2 can be set to 1 (e.g., Burger and Lynch 1995). The key 
theoretical tool for studying phenotypic evolution is the 
Lande equation (Lande 1976a), whose univariate version 
reads 

Az f = a\p t , (1) 

where Az t is the change in mean phenotype after one gen- 
eration of selection, and /? ( = d{\nw t ) / dz t is the selection 
gradient at time t, that is, the derivative of log mean fitness 
w t with respect to the mean phenotype. Note that eqn (1) 
is analogous to the univariate breeder's equation 
Az t = h 2 S t , where S t = o 2 fi t = cov(w t ,z f ) is the selec- 
tion differential. A rate of change in haldanes can be 
obtained by standardizing with a p , yielding 

^ = (2) 

where = cov(w f , z t )/a p is the variance -standardized 
selection gradient (Lande and Arnold 1983; Hereford et al. 
2004). 

For multiple traits, the structure of phenotypic variation 
is summarized by the matrix P, whose diagonal entries con- 
tain the phenotypic variances of the individual traits, and 
whose off-diagonal entries contain the phenotypic covari- 
ances. In the standard model, P = G+E, where G is the 
(additive) genetic covariance matrix and E the matrix of 
environmental variances and covariances. The multivariate 
version of Lande's equation is 

Ait = Gp t , (3) 

where, for n traits, z t = (z l t , . . .,z„.,)' is the vector of 
mean trait values (with ' denoting transposition) and 
Pt = (dw t /dzij, ■ . .jdwt/dznj)' is the multivariate selec- 
tion gradient, which points in the direction of steepest 
ascent on the fitness landscape. The response to selection is 
also influenced by the structure of genetic variation speci- 
fied in the G-matrix. In particular, genetic correlations can 
cause the response to selection to show a bias toward trait 



combinations with high genetic variation (see Fig. 2 below; 
for an introduction to the geometric aspects of multivariate 
selection, see Walsh and Blows 2009). 

The structure of multivariate genetic variation is often 
analyzed in terms of the eigenvectors of the G-matrix (as in 
a principal component analysis). The eigenvectors (princi- 
pal components) can be viewed as composite traits (linear 
combinations of the original traits) that are genetically un- 
correlated (i.e., their covariances are zero) and whose 
genetic variances are given by the corresponding eigen- 
values. Graphically, if the distribution of breeding values 
(i.e., the average contribution of an individual to the phe- 




Traiti 




Trait i 



Figure 2 Illustration of adaptation involving two genetically correlated 
traits. (A) Adaptation after a sudden environmental change; the new 
optimum 0, is constant. Gray lines illustrate the fitness surface, defined 
by the matrix w. The distribution of breeding values defined by the G- 
matrix is illustrated by the black ellipse, whose center is the mean phe- 
notype z t and whose axis is the eigenvectors of G. The initial response 
to selection is biased toward the leading eigenvector, that is, the 
genetic line of least resistance (Schluter 1996). (B) Adaptation to a mov- 
ing optimum. Gray circles show the fitness landscape at four different 
points in time. Black ellipses show the corresponding positions of the 
population (represented by the G-matrix). The insets at the top show 
the leading eigenvector of G, k\, the selection gradient /f t and the 
response to selection Az t at time points 1, 2, and 3, respectively. 
Because the initial response is biased toward the leading eigenvector, 
the population 'rises' above the line of the moving optimum (i.e., the 
flying-kite effect; Jones et al. 2004). This rise comes to a halt as the ten- 
dency to follow the line of least resistance is balanced by the selection 
gradient, resulting in horizontal movement of the population. 
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notype of its offspring) is multivariate Gaussian, isoclines 
of this distribution can be represented by ellipses (or 
higher-dimensional ellipsoids), with axes given by the ei- 
genvectors and their lengths proportional to the roots of 
the eigenvalues (Fig. 2). The major axis of such an ellipse 
(i.e., the leading eigenvector of the G-matrix) represents 
the trait combination with a maximum of genetic variation. 
Is has been called g max or the genetic line of least resistance 
(Schluter 1996). Eigenvectors with small (or zero) eigen- 
values represent trait combinations with little (or no) 
genetic variation, into which evolution is severely con- 
strained (Hansen and Houle 2008; Gomulkiewicz and Ho- 
ule 2009; Kirkpatrick 2009; Walsh and Blows 2009; Chevin 
2013). More generally, it is also possible to calculate the 
amount of variation along any direction of the phenotypic 
space (Hansen and Houle 2008; Gomulkiewicz and Houle 
2009). For the pros and cons of multivariate analysis in 
quantitative genetics, see Houle et al. (2002), Mezey and 
Houle (2003), Pigliucci and Kaplan (2006), Blows (2007), 
Walsh and Blows (2009), Berner (2012) and the commen- 
taries to Blows (2007) in volume 20:1 of the Journal of Evo- 
lutionary Biology. 

Phenotypic plasticity 

Phenotypic plasticity in quantitative traits is usually charac- 
terized by reaction norms, which give the phenotype as a 
function of an environmental variable. When different 
genotypes have different reaction norms, plasticity is itself 
evolvable. While the evolution of plasticity can be modeled 
in different ways (Via and Lande 1985; De Jong 1995, see 
also Box 1 in Chevin et al. 2013), most of the models 
reviewed here focus on linear reaction norms and treat 
their slope and elevation as quantitative traits (e.g., Lande 
2009). The majority of models have studied plasticity in 
single traits only (but see Gavrilets and Scheiner 1993; 
Draghi and Whitlock 2012), even though the G-matrix is 
known to be sensitive to environmental conditions (e.g., 
Tonsor and Scheiner 2007; Husby et al. 2011). 

Population dynamics 

Models of evolutionary rescue assume that the intrinsic 
population growth rate depends on the degree of adapta- 
tion, that is, on mean absolute fitness. Regardless of poten- 
tial density dependence, a population will decline if the 
average number of offspring per individual drops below 1. 
That is, eventually, population size N is likely to follow 

N t+ i = w t N t . (4) 

As shown in Appendix 1, the mean fitness w t is generally 
reduced by two kinds of genetic load (Lande and Shannon 
1996; Chevin 2013): a standing load due to phenotypic vari- 
ation and a lag load (Maynard and Smith 1976) due to 
deviations of the mean phenotype from the optimum (also 



called selection load). In many models, survival or extinc- 
tion of the population depends primarily on the lag load. A 
crucial point is that population dynamics depend on the 
mean fitness (eqn 4), whereas evolutionary change depends 
on the fitness gradient (eqns 1 or 3). Another way of saying 
this is that population dynamics depend on absolute fitness 
and evolution on relative fitness (Bell 2013). The relation- 
ship between these two quantities is determined by the fit- 
ness function: a given fitness gradient can be associated 
with a higher mean fitness under strong selection than 
under weak selection (Fig. IB). This point will be essential 
in our discussion of sustainable evolutionary rates (see 
below). 

Some predictions from the various models reviewed in 
this study are summarized in Table 1 . 

Adaptation of a single quantitative trait 

Sudden environmental change 

In the sudden-change scenario, a population that is well 
adapted to its environment is displaced from the fitness 
peak by a sudden shift of the optimum. Phenotypic evolu- 
tion is relatively straightforward: the mean phenotype will 
approach the new optimum exponentially (because the fit- 
ness gradient decreases in the vicinity of the optimum) 
(Lande 1976b). The key question is whether evolution is 
fast enough in cases where, immediately after the environ- 
mental change, the population mean fitness is less than 1. 
In this case, the population size will initially decline, setting 
off a 'race' between adaptation and extinction. Gom- 
ulkiewicz and Holt (1995) showed that evolutionary rescue 
is possible only if the initial maladaptation after the envi- 
ronmental change is not too large and the initial popula- 
tion size is high. 

Gradual environmental change 

The situation is quite different if the optimum changes 
gradually rather than suddenly. In the simplest case, the 
optimum increases linearly at rate k. This model has been 
analyzed by Lynch et al. (1991), Lynch and Lande (1993), 
and Burger and Lynch (1995) and later been extended by 
various authors (see below). An excellent summary is given 
in the study by Burger and Lynch (1997). As the behavior 
of this model is highly instructive, we will describe it in 
some detail (see also Appendix 1). 

Assume again that the original population is well 
adapted. As the optimum starts moving, selection becomes 
gradually stronger (see eqn A4). Consequently, the popula- 
tion will initially evolve slowly, and the lag between the 
optimum and the population mean phenotype will increase 
(the population 'slips off the fitness peak). However, as the 
distance to the optimum increases, so does the selection 
gradient, until finally a state of dynamic equilibrium is 
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Table 1 . A summary of theoretical predictions for models of adaptation to environmental change 



Mode of 

Env. 

Change 



Phenotypic Evolution 



Survival/Extinction 



Effects of Plasticity 



Sudden Change 

Single Trait Change in mean phenotype described by 
univariate Lande equation (eqn 1). 
Exponential approach to new optimum 
(Gomulkiewicz and Holt 1995). 



Multiple Traits 



Gradual Change 



Change in mean phenotype described by 
multivariate Lande equation (eqn 3). 
Trajectory to new optimum biased 
toward genetic line of least resistance 
(Schluter 1996). Lag load decreases 
roughly exponentially (Chevin 2013). 

Population follows the optimum with a 
constant lag (Lynch and Lande 1993). 
Trait correlations can induce permanent 
maladaptation in traits under stabilizing 
selection ('flying-kite effect'; Jones et al. 
2004). Genetic variance increases 
(Burger and Lynch 1995; Burger 1999). 



Random Change 

Population's ability to track the optimum 
increases with autocorrelation of 
fluctuations (Lande and Shannon 1 996; 
Chevin 2013). Autocorrelated 
fluctuations increase genetic variance, 
whereas uncorrelated fluctuations do 
not (Burger 1999). 
Spatial Heterogeneity 

Spatial heterogeneity constrains 
adaptation. Trait interactions can induce 
'counter-gradient' dines, causing traits 
to evolve away from the optimum 
(Duputieet al. 2012). 



Maximal amount of environmental change 
the population can handle depends on 
width of fitness function, intrinsic growth 
rate, initial population size, and genetic 
variation (Gomulkiewicz and Holt 1995) 

As for single trait. Extinction risk depends on 
genetic variance in direction of selection 
(Gomulkiewicz and Houle 2009). 



Critical rate of environmental change (eqn A6) 
increases with genetic variance in direction 
of moving optimum and with intrinsic 
growth rate and is maximal at intermediate 
strength of selection (Lynch and Lande 
1993; Burger and Lynch 1995; 
Gomulkiewicz and Houle 2009), see Fig. A1 . 



Extinction risk elevated if fluctuations are 
uncorrelated and occur in directions with 
strong selection and high genetic variation 
(Chevin 2013). 



Population growth is maximized at 
intermediate dispersal rates. Critical rate of 
environmental change increases if spatial 
selection gradient is aligned with direction 
of abundant genetic variation and weak 
stabilizing selection (Duputie et al. 2012). 



Approach to new optimum is facilitated 
by temporary increase in phenotypic 
plasticity and a concomitant release of 
hidden genetic variation (Lande 
2009). Plasticity reduces extinction risk 
(Chevin and Lande 2010). 

No models available for the context of 
environmental change. 



Adaptive plasticity reduces the 
perceived speed of environmental 
change. => increases critical rate of 
change, decreases phenotypic lag, 
and decreases rate of genetic 
evolution. Effects may be 
counteracted by costs of plasticity 
(Chevin et al. 2010). 

Strong plasticity increases extinction 
risk if environmental cues are 
unreliable (Reed et al. 2010). 
Predictable fluctuations can select for 
increased plasticity. 



If plasticity is expressed before (after) 
migration, it increases (reduces) 
migration load and can decreases 
(increase) species ranges. Expressed 
plasticity increases near range limits 
(Chevin and Lande 2011; Thibert- 
Planteand Hendry 2011). 



reached, at which the rate of evolution exactly matches the 
rate of environmental change (see Fig. 1A and eqn A5). 
Whether or not the population survives depends on the 
mean fitness at this distance from the optimum (i.e., on the 
lag load, which is approximately proportional to k 2 ; Lande 
and Shannon 1996). One can thus calculate a critical rate of 
environmental change k CTit (eqn A6), which is the maximal 
rate of change the population can handle. If the environ- 
ment changes faster than /c crit , the lag load becomes so large 
that the population can no longer maintain itself. Extinc- 
tion usually follows quickly, because the reduction in pop- 
ulation size leads to a loss of genetic variation, which 
further undermines the population's ability to adapt. 



Thus, in contrast to the sudden-change scenario, evolu- 
tionary rescue in a gradually changing environment 
requires that the population maintain a positive growth 
rate at all times. This is a consequence of the 'relentless' 
movement of the optimum, which means that a population 
that has fallen behind in the race will get no chance to catch 
up. It also is noteworthy that extinction in this model usu- 
ally is not due to a lack of genetic variance (except in the 
final phases of the collapse), nor due the classical 'cost of 
selection' (i.e., the required number of selective deaths, 
Haldane 1957). Rather, the population dies out because all 
individuals (not just the less adapted ones) have low fit- 
ness. 
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The critical rate of environmental change is directly pro- 
portional to the additive genetic variance and the square 
root of the maximal population growth rate (see eqn A6). 
The dependence on the width of the fitness landscape - or 
conversely, the strength of stabilizing selection - is more 
complex: as shown in the first row of Fig. Al, for constant 
o^> fccrft 1S maximal at small to intermediate values of the 
parameter 'V/ which measures the effective width of the 
fitness function. In other words, the population can sup- 
port the fastest environmental change if stabilizing selec- 
tion is strong, but not too strong. The drop-off in k clit at 
low or high values of V s can be explained by the two kinds 
of genetic load introduced above. On the one hand, very 
strong selection (i.e., in a steep and narrow fitness land- 
scape; small Vj) induces a high standing load, which 
reduces the realized growth rate and diminishes the ability 
of the population to tolerate environmental change. On the 
other hand, sufficiently weak selection in combination with 
a moving optimum increases lag load, because the popula- 
tion will follow the optimum at a greater distance. This 
somewhat counter-intuitive result is due to the fact that, 
on a flatter fitness landscape, reaching a given selection gra- 
dient requires a larger decrease in mean population fitness 
(see above and Fig. 1). In other words, whereas strong 
selection keeps the population close to the optimum at 
high mean fitness, weak selection, precisely because it is 
ineffective, allows the population to slip farther off the fit- 
ness peak. Therefore, weak selection in combination with a 
constantly moving optimum represents a 'slippery slope' 
that can be very dangerous for population survival (see dis- 
cussion in Burger and Lynch 1995 and Huey and Kingsolv- 
er 1993). Biirger and Lynch (1995) also showed that the 
critical rate of change is further decreased by genetic drift 
of the mean phenotype in small populations and by sto- 
chastic fluctuations of the optimum around the linear 
trend (see also Bjorklund et al. 2009). 

In many quantitative-genetic models, the additive 
genetic variance a 2 is assumed to be constant. Over short 
timescales, this may be approximately true, but over longer 
timescales, a 2 is itself subject to evolutionary change, and 
it is this fact that makes expressions for k clit (such as 
eqn A6) 'deceptively simple' (Burger and Lynch 1995). 
Explaining the evolution and maintenance of genetic varia- 
tion is one of the perennial problems in theoretical popula- 
tion genetics, and no fully satisfactory model has as of yet 
been found (Barton and Turelli 1989; Burger 2000; Barton 
and Keightly 2002; Johnson and Barton 2005; Hill 2010). 
Before the environmental change, the population may be 
assumed to be at mutation-selection-drift balance, for 
which several approximations have been developed (Lande 
1976a; Turelli 1984; Burger 2000; Alvarez-Castro et al. 
2009). In the second row of Fig. Al, we follow Burger and 
Lynch (1995) by showing the predicted values of k ait (in 



units of the phenotypic standard deviation a p , see below) 
when a 2 is chosen according to the so-called stochastic 
house-of-cards approximation (Burger et al. 1989). Doing 
so takes into account that populations under weak selec- 
tion have higher genetic variance, which may offset the 
negative effects of weak selection on the lag load (see 
above) and lead to a positive relationship between the 
width of the fitness landscape and k cA (see Huey and King- 
solver 1993). However, this is still not the whole story, 
because once the optimum starts moving, a 2 is expected to 
increase. This increase is mainly due to the rise in fre- 
quency of previously rare alleles, and it is strongest in large 
populations (Burger 1999): for example, under standard 
values of mutational and selectional parameters, cr 2 
increases up to 4-fold in populations with N e > 5000. In 
contrast, selection has little impact on a 2 if N e < 200—300 
(Burger 1999), which might explain why genetic variances 
usually do not increase in artificial selection experiments, 
as noted by (Johnson and Barton 2005). A useful upper 
limit for the genetic variance in small populations 
(N e < 500, Burger and Lynch 1995) is the neutral expecta- 
tion 2 V m N c , where V m is the input of genetic variance from 
new mutations (a typical value is V m = O.OOler;, Lande 
1976a; Lynch 1988). In summary, evolution of the genetic 
variance may increase the prospects of population survival, 
but mostly in large populations. It should be noted, 
though, that the increase in variance takes time and may 
come too late for populations subject to strong environ- 
mental change. 

Fluctuating selection 

In addition to sudden or gradual changes, most environ- 
ments are subject to stochastic fluctuations. We have 
already seen that superimposing fluctuations on a linear 
trend in the optimal phenotype increases population 
extinction risk and decreases the critical rate of environ- 
mental change A: crit (Burger and Lynch 1995). Here, we 
briefly discuss the effects of fluctuations around a constant 
mean. Uncorrelated fluctuations (white noise) in the opti- 
mal phenotype resemble a sudden-change scenario that is 
repeated each generation. Such fluctuations can incur 
strong selection, but the responses of the population wiE 
not add up to large changes over longer timescales (Ginge- 
rich 1983; Gibbs and Grant 2006). In addition, genetic 
responses to selection in one generation are likely to be 
maladaptive in the next generation, and therefore, the lag 
load will be high (Lande and Shannon 1996; Burger 1999; 
Chevin 2013). Consequently, uncorrelated fluctuations do 
not lead to a significant increase in genetic variance relative 
to constant stabilizing selection (Burger 1999). An excep- 
tion exists, however, if a species possesses dormant stages 
such as seeds or resting eggs or if generations are overlap- 
ping but selection acts only on juveniles. In these cases, the 
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'storage effect' allows the maintenance of genetic polymor- 
phism and, hence, high levels of variation (Chesson and 
Warner 1981; Hairston et al. 1996). Environmental fluctu- 
ations can also select for phenotypic plasticity, provided 
the state of the environment can be assessed by a reliable 
cue (Tufto 2000), or for bet-hedging, if there is no such cue 
(Svardal et al. 2011). 

In contrast to uncorrected fluctuations, autocorrelated 
fluctuations are more similar to the gradual-change 
scenario, and a population with sufficient genetic variance 
can follow the optimum and maintain high fitness 
(Charlesworth 1993; Lande and Shannon 1996; Chevin 
2013). Consequently, autocorrelated fluctuations can lead 
to significant increases in genetic variation (Burger 1999). 

Adaptation of multiple correlated traits 

When several traits are under selection, the above analyses 
need to be extended to account for the effects of genetic 
correlations. As mentioned above, genetic correlations tend 
to bias the phenotypic response to selection toward the 
leading eigenvector of the G-matrix, g max (the 'genetic line 
of least resistance'; Schluter 1996). In the sudden-change 
scenario, an evolving population will still reach the new 
optimum, although not along the most direct path 
(Fig. 2A). While the optimum is approached, the lag load 
decreases as a sum of exponential terms, with rates given by 
the eigenvalues of the matrix of selection responses (Chevin 
2013). Adaptation is fastest and evolutionary rescue is most 
likely if the angle between the direction of selection and 
gmax is small (Gomulkiewicz and Holt 1995; Schluter 
1996). 

Under gradual environmental change, selection for a 
moving optimum may cause permanent maladaptation 
of traits (or trait combinations) that are under pure sta- 
bilizing selection (i.e., orthogonal to the direction of the 
optimum). As illustrated in Fig. 2B, the initial response 
to selection is biased toward g max , causing the population 
to rise above the line of the moving optimum, a phe- 
nomenon that has been termed the 'flying-kite effect' 
(Jones et al. 2004). Eventually, the rise comes to a halt, 
as stabilizing selection in the respective direction 
increases, and the population's trajectory continues in 
parallel to that of the optimum. Again, population sur- 
vival will depend on the lag load at this steady state. A 
critical rate of environmental change can be calculated in 
analogy to the univariate case (see Appendix 2). It 
depends not only on the shape of the fitness landscape, 
but also on the direction of the optimum and the struc- 
ture of the G-matrix. In particular, the critical rate is 
high if the optimum moves in parallel to g max , and it is 
lowest if the optimum moves in a direction of low 
genetic variation (see Hellmann and Pineda-Krch 2007 



for graphical illustrations and a discussion of the conse- 
quences for conservation biology). 

As in the univariate case, many studies assume that the 
G-matrix is roughly constant over the timescale of inter- 
est. Evolution of the G-matrix has been studied in a 
recent series of papers by Jones, Arnold and Biirger 
(Jones et al. 2003, 2004 2007, 2012; for review see Arnold 
et al. 2008). In accordance with previous studies (Barton 
and Turelli 1987; Burger and Lynch 1995; Jones et al. 
2004), Jones et al. (2012) found that, irrespective of the 
mode of environmental change (gradual, episodic, sto- 
chastic), genetic variance increases in the direction of 
environmental change. While this facilitates the response 
to selection, the phenotypic lag also induces a skew in 
the distribution of breeding values (unfit phenotypes 
'trailing behind'), which restrains the response to selec- 
tion. Generally, the two phenomena do not offset each 
other (Jones et al. 2012), requiring inspection for every 
individual case. These results highlight the need for cau- 
tion when iterating the Lande equation or interpreting 
G's eigenvalues (Kirkpatrick 2009). Under pure stabilizing 
selection, the G-matrix tends to align itself with the fit- 
ness landscape, that is, genetic variance is highest in 
directions with weak selection. G depends, however, also 
on the distribution of new mutations, that is, the M- 
matrix (Jones, et al. 2003, 2007), and on gene flow (Guil- 
laume and Whitlock 2007; Franks et al. 2014). 

Genetic basis of adaptation 

The quantitative-genetic models we have considered so far 
are most accurate if adaptation is based on a large number 
of loci with small individual effects. In this section, we 
briefly discuss several issues that arise when this assump- 
tion is relaxed. 

The first question is how the rate of adaptation is 
affected by alleles of large effect. If the same total progress 
toward the optimum can be made by the fixation of either 
a single allele of large effect or many alleles with small 
effects, adaptation will be faster in the former case, because 
selection on the large alleles is more effective (Gom- 
ulkiewicz et al. 2010; for the same result in a different con- 
text, see also Gavrilets et al. 2007; Rettelbach et al. 2011). 
In Appendix 3, we calculate the rate of phenotypic evolu- 
tion due to the fixation of a major allele and show that it 
can be quite high, at least while the allele is at intermediate 
frequency. For quantitative traits that are determined by a 
combination of small- and large-effect loci, Gomulkiewicz 
et al. (2010) showed that adaptation is fastest when both 
classes of loci are evolving. For the same situation, Chevin 
and Hospital (2008) demonstrated that 'background'-adap- 
tation from minor loci, by successively reducing the selec- 
tive advantage of a large-effect allele, can significantly affect 
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its trajectory, and even prevent fixation. The exact outcome 
crucially depends on the initial allele frequency, the dis- 
tance from the optimum, and the amount of genetic varia- 
tion provided by the minor loci. 

Another question is, however, how likely beneficial alleles 
with large effect are in the first place. In a multivariate con- 
text, Fisher (1930) used his classical 'geometric model' to 
argue that alleles (i.e., mutations) with large effect that 
pleiotropically affect multiple traits are most likely to be 
deleterious. As pointed out by Kimura (1983), however, 
Fisher neglected the fact that, among beneficial mutations, 
the few mutations with large effect have a higher fixation 
probability than the more common mutations with small 
effects. In the last two decades, numerous theoretical stud- 
ies have developed predictions for the distribution of phe- 
notypic and fitness effects of both new and fixed mutations 
(e.g., Martin and Lenormand 2006a, 2008; Keightley and 
Eyre-Walker 2007; Yeaman and Whitlock 2011), and many 
models have concluded that the role of mutations with 
major effects in adaptation is surprisingly large (reviewed 
by Orr 2005). However, almost all of these models have 
considered a sudden-change scenario. Under gradual envi- 
ronmental change, results might be very different. In par- 
ticular, Collins et al. (2007) and Kopp and Hermisson 
(2007, 2009a,b) showed that a slowly moving optimum 
favors adaptation by small mutations. 

Finally, many authors have studied adaptation and evo- 
lutionary rescue from a single large mutation. As these 
models usually do not refer to quantitative traits, we only 
point out the relevant literature: for the probability of evo- 
lutionary rescue, see Gomulkiewicz and Holt (1995); Holt 
and Gomulkiewicz (1997); Orr and Unckless (2008); Uec- 
ker and Hermisson (2011); Martin et al. (2013); for the fix- 
ation probability of a new mutation in a changing 
environment, see Uecker and Hermisson (2011); Kirkpa- 
trick and Peischl (2013); Martin et al. (2013); and for the 
probability of adaptation from standing genetic variation 
versus new mutations, see Hermisson and Pennings (2005); 
Martin et al. (2013). 

Maximal sustainable rates of evolution? 

A well-known prediction from the models by Lynch and 
Lande (1993) and Burger and Lynch (1995) is that of a 
'maximal sustainable rate of evolutionary change' on the 
order of 0.1 haldanes or less. This value is simply a ballpark 
estimate of the critical rate of gradual environmental 
change, A: crit (eqn A6), scaled by the phenotypic standard 
deviation and parametrized with realistic parameter values 
(see Appendix 1 and Fig. Al). As, at the dynamic equilib- 
rium, the population follows the optimum with a constant 
lag, the rates of environmental and phenotypic change are 
'formally equivalent' (Burger and Lynch 1995). For clarity, 



we will denote the rate of phenotypic change in haldanes by 
Kcrit = KrulOp (eqn A7). 

Barrett and Hendry (2012) note that it is 'tempting' to 
use K crit = 0.1 as a benchmark for empirically observed 
evolutionary rates, the idea being that rates near or above 
this value might be cause for concern because they are not 
'sustainable' (see also Hendry and Kinnison 1999). Based 
on earlier meta-analysis (Hendry and Kinnison 1999; Kinn- 
ison and Hendry 2001; Hendry et al. 2008), these authors 
conclude that most rates of change are below 0.1. In con- 
trast, Gingerich (2009) argued that evolutionary rates on 
the order of 0.1 and 0.3 haldanes are common, but his 
analysis relied on an interpolation technique (log-rate-log- 
interval plots) that is sensitive to measurement error when 
real rates of change are small (Hunt 2012 and below). Bar- 
rett and Hendry (2012) also warn, however, that theory- 
derived critical rates rely on 'many unrealistic assumptions, 
such as perpetual persistence under constant environmen- 
tal change' and that 'critical rates for natural populations 
over time frames of conservation interest could be very dif- 
ferent'. 

There are several points to be made here (see also 
Appendix 4). First, and obviously, a universal K clit of 0.1 
haldanes cannot be more than a rule of thumb. Critical 
rates may be higher under strong selection, high heritabili- 
ties and in large populations (Fig. Al). Second, some of the 
reasons for population extinction found by Burger and 
Lynch (1995) - such as random but autocorrelated fluctua- 
tions in genetic variance - are, indeed, mainly a long-term 
concern under sustained environmental change. Third, 
however, the critical rate in eqns (A6) and (A7) is simply 
equivalent to the (instantaneous) rate of evolutionary 
change that can be achieved without a decrease in popula- 
tion size, as a function of (i) the genetic variance, (ii) the 
reproductive capacity of the population, and (iii) the shape 
of the fitness landscape (see discussion of mean fitness vs. 
fitness gradient above). Faster evolution is possible tempo- 
rarily, but only at the cost of a reduction in population size. 
To quantify this effect, in Appendix 4, we estimate maximal 
rates of environmental and phenotypic change when allow- 
ing modest population decline over a limited time frame 
(e.g., the population is to maintain a minimal size of 50 
individuals for 50 generations). As shown in Fig. A2, this 
provision leads to modest increases in K crit in large popula- 
tions (typically around 30%), whereas the effect in small 
populations is negligible (in particular, in the light of the 
stochastic variations discussed below). In summary, K crit is, 
indeed, likely to often be around or below 0.1 haldanes. 
Faster observed rates may be a sign that the population is 
under stress (e.g., the well-known example of Darwin's 
finches during a drought, where beak-size increased by 0.66 
standard deviations, but 85% of the population died; Grant 
and Grant 2006) or may indicate that part of phenotypic 
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change is due to plasticity (see below). Temporarily high 
rates of change may also be achieved by the fixation of a 
large-effect mutation (Appendix 3). 

Small maximal rates of phenotypic change also raise sta- 
tistical issues (Appendix 5; see also Hendry and Kinnison 
1999): detecting a difference of 0.1 standard deviations 
between two populations requires very large sample sizes 
(e.g., almost 800 per population for 50% power in a two- 
sample f-test with « = 0.05). On the other hand, differences 
of this magnitude can easily be created by sampling effects 
(Figs 5, 6; Kinnison and Hendry 2001; Hunt 2012). Indeed, 
the mean absolute differences in units of phenotypic stan- 
dard deviations between two samples of size n drawn from 
the same population is 2/ ' y/nk (Hunt 2012), which equals 
0.113 for n = 100. In finite populations, similar effects 
occur due to genetic drift and environmental variance 
(even if the whole population is sampled). The variance of 
the mean phenotype due to genetic drift is a 2 , /N e , with N e 
being the effective population size (Lande 1976a). By a cal- 
culation analogous to the one in Hunt (2012), the mean 
generation-to-generation rate in haldanes due to drift is 
2^h 2 /(nN e ), which is 0.025 for N e = 1000 and h 2 = 0.5. 
Similarly, the contribution of environmental variance (i.e., 
genotype-independent random variation in individual phe- 
notypes) to the mean rate of phenotypic change is 
2^(1 - h 2 )/(nN) (with N being the census population 
size). Together, these two sources of variation may domi- 
nate the generation-to-generation changes in the mean 
phenotype of small populations (Appendix 5, Figs A6, A7). 
In summary, maximal sustainable rates of evolutionary 
change might often be of the same order than various 
sources of stochastic noise, something which should be kept 
in mind when interpreting evolutionary rates measured 
over short timescales. 

The role of phenotypic plasticity 

So far, we have only considered genetic adaptation. How- 
ever, many observed responses to climate change are likely 
to be plastic (Gienapp et al. 2008; Hendry et al. 2008; Mer- 
ila 2012), and assessing the relative importance of plastic 
and genetic changes is precisely the aim of this special issue 
of Evolutionary Applications. Yet, in its basic form, the 
question is empirical and cannot be answered by theory 
alone. While quantitative genetic models can make some 
tentative predictions about the maximal rates of genetically 
based evolution (see above), it seems impossible to make 
general statements about the range and scope of plasticity. 
Here, we will instead focus on reviewing models that inves- 
tigate the interaction between plasticity, population 
dynamics, and genetic evolution. Because several important 
aspects have already been reviewed elsewhere (Ghalambor 
et al. 2007; Chevin et al. 2013), our treatment can be short. 



Ecological models have investigated the effect of plastic- 
ity on population stability and extinction risk in the 
absence of evolution. Community models including 
so-called trait-mediated indirect effects (Werner and Pea- 
cor 2003) frequently find that phenotypic plasticity 
mediated by species interactions (e.g., inducible defenses 
against predators; Tollrian and Harvell 1999) can stabilize 
population dynamics, even though such a stabilizing influ- 
ence is not universal (Kopp and Gabriel 2006). If plasticity 
increases the range of conditions under which a commu- 
nity is stable, it reduces the risk of species extinctions after 
an arbitrary environmental change ('plastic rescue'; Kov- 
ach-Orr and Fussmann 2013). 

For a single population, Reed et al. (2010) studied the 
impact of phenotypic plasticity on population extinction 
risk in a randomly fluctuating environment. They found 
that adaptive plasticity decreases extinction risk, unless the 
magnitude of plastic responses exceeds an optimal level set 
by cue reliability (strong responses to unreliable cues tend 
to be harmful). Chevin et al. (2010) included phenotypic 
plasticity into the moving-optimum model of Lynch and 
Lande (1993) and Burger and Lynch (1995). Assuming a 
linear reaction norm with slope less than one, plasticity 
essentially reduces the speed of environmental change per- 
ceived by the population. Plasticity thus increases the criti- 
cal rate of environmental change k ait that separates 
population survival from extinction. In consequence, it 
increases the maximal rate of phenotypic change, while 
simultaneously decreasing the rate of genetic evolution. 
This effect may be reversed at high levels of plasticity if 
plasticity itself is costly (and hence, reduces the mean fit- 
ness of the population). 

Gienapp et al. (2013) recently applied both the Chevin 
et al. (2010) and the Burger and Lynch (1995) model to 
anticipate evolution of egg-laying dates in great tits from a 
well-studied Dutch population. Egg-laying date in this spe- 
cies is a phenotypically plastic trait that depends on spring 
temperature and is selected to coincide with the peak in 
caterpillar abundance. Using various modeling techniques, 
the authors show that, despite plasticity, global warming 
will create a mismatch between the optimal and realized 
egg-laying dates, which might threaten population persis- 
tence unless it can be closed by genetic evolution. By focus- 
ing on the predicted mismatch, the authors were able to 
parametrize the Burger and Lynch (1995) model (i.e., 
eqn A6), even though this model was not built to deal with 
plasticity. They conclude that, even under a mild climate- 
change scenario, the predicted rate of environmental 
change (from the point of view of the population) is close 
to the theoretical maximal sustainable rate. To parametrize 
the Chevin et al. (2010) model, Gienapp et al. (2013) 
assumed that both optimal and realized egg-laying dates 
correlate with mean spring temperature (measured between 
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mid-March and mid-April). Although the Chevin et al. 
(2010) model seems to be more suitable for the analysis of 
a plastic trait, its results are less plausible than those 
obtained from the Burger and Lynch (1995) model. In par- 
ticular, the model predicts that population survival will be 
facilitated by fast environmental change. The authors argue 
that this counterintuitive prediction is an artifact, which 
arises because, with faster temperature increase, mean 
spring temperature becomes less and less correlated with 
the true causal variable determining optimal egg-laying 
date. This highlights the general problem that, frequently, 
the variables we can measure are just proxys for one or 
more causal factors. If the proxy is bad, any model will per- 
form poorly. Despite these issues, the study by Gienapp 
et al. (2013) is exemplary in its combined use of long-term 
empirical data, climate-change predictions, and models for 
future optimal and realized behavior. 

In the following, we briefly review models in which plas- 
ticity can itself evolve. Conditions for the evolution of plas- 
ticity are fairly well understood. Plasticity is adaptive if 
individuals encounter different environmental conditions 
that favor different phenotypes and that can be assessed by 
a reliable cue (e.g., Tollrian and Harvell 1999; Ghalambor 
et al. 2007). Its evolution may be limited by functional 
constraints, unreliable cues (Tufto 2000) and costs for the 
necessary sensory and developmental machinery (DeWitt 
et al. 1998; van Buskirk and Steiner 2009). More recently, 
however, phenotypic plasticity has been advocated as not 
only a product, but also a driver of genetic evolution 
(West-Eberhard 2003; for recent reviews, see Ghalambor 
et al. 2007; Pfennig et al. 2010; Wennersten and Forsman 
2012; Wund 2012). The basic idea is that new phenotypes 
first appear as a result of environmental induction and only 
later are fixed via 'genetic assimilation' or 'genetic accom- 
modation'. Here, genetic assimilation corresponds to a loss 
of plasticity, such that expression of the phenotype 
becomes independent of environment cues. Genetic 
accommodation is a more general 'fine-tuning' of the novel 
phenotype via changes in allele frequencies, potentially 
facilitated by a release of hidden genetic variation (Hermis- 
son and Wagner 2004; Moczek 2007; for more conceptual 
discussion, see West-Eberhard 2005; Crispo 2007; Ghalam- 
bor et al. 2007). The more ambitious versions of this 
hypothesis - that environmental induction can be at the 
basis of 'evolutionary novelties' (West-Eberhard 2003; Pig- 
liucci et al. 2006; Uller and Helantera 2011) - appears 
unaccessible to classical population-genetics modeling. 
Here, we focus on the less far-reaching question of the 
role of plasticity in the evolution of existing quantitative 
traits. 

Phenotypic plasticity has traditionally been viewed as 
delaying genetic evolution. This is certainly true if plasticity 
is sufficient to ensure continued high fitness of a popula- 



tion in a changing environment. However, there are other 
scenarios in which plasticity may, indeed, speed up or facil- 
itate genetic change. A simple case is the Baldwin effect 
(Baldwin 1896; Crispo 2007), where plasticity (specifically, 
learning) allows a population to survive in a new or chan- 
ged environment, thereby enabling future genetic adapta- 
tion (for models, see Ancel 1999; Pal and Miklos 1999; 
Ancel 2000; Paenke et al. 2007). Furthermore, plasticity 
can influence the course of evolution by bringing a popula- 
tion into the domain of attraction of a specific adaptive 
peak. The probability of a peak shift is highest if plasticity 
is of intermediate strength (Price et al. 2003). Both mecha- 
nisms may play a role in biological invasions as well as 
adaptation to climate change. 

Recently, Lande (2009) proposed a simple model for the 
role of plasticity in adaptation to an abrupt environmental 
shift. He considered the evolution of a quantitative trait 
that is determined by linear reaction norms. That is, for 
each individual, the trait value is a linear function of an 
environmental variable, with genetic variation in the slope 
and intercept of this function (see also Gavrilets and Schei- 
ner 1993). Under the original conditions, a modest level of 
plasticity (i.e., an intermediate slope of the reaction norm) 
is favored in a slightly fluctuating environment with con- 
stant mean and imperfect cues. At this stage, reaction-norm 
slope varies between individuals, but the mean phenotype 
is relatively homogeneous (canalization). When the mean 
environment changes, genetic variance is increased due to 
differential plastic responses (decanalization), and selection 
favors individuals with steep reaction norms, which can 
best adjust to the new conditions. That is, the population 
evolves toward the new optimum via the evolution of 
increased plasticity, allowing high rates of phenotypic 
change. Subsequently, the reaction norm intercepts 
increase and slopes decrease, again reaching the optimal 
degree of plasticity in the new environment (genetic assimi- 
lation). Chevin and Lande (2010) added population 
dynamics to this model and showed that evolving plasticity 
strongly increases the probability of evolutionary rescue 
after a sudden environmental change. 

What determines adaptive potential? 

Ideally, we would like to be able to predict which species 
have the potential to adapt to rapid climate change (Wil- 
liams et al. 2008; Huey et al. 2012). Obviously, phenotypic 
plasticity will help (see above), but theory can say little 
more than that. With regard to genetic adaptation, the 
adaptive potential depends most directly on the genetic 
variation that is available in the direction of selection. In 
addition, we may also ask what kind of genetic architec- 
tures and evolutionary histories facilitate rapid adaptation. 
We will discuss these two issues in turn. 
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Genetic variance and genetic constraints 

For single traits, a short-term measure of adaptive potential 
is given by the additive genetic variance (see eqn 1), and a 
lack of such variance corresponds to a genetic constraint 
(i.e., adaptive potential and genetic constraints are two 
sides of the same coin). An absolute constraint is present if 
genetic variance is zero, and a relative constraint if it is low. 
Gomulkiewicz and Houle (2009) pointed out that if adap- 
tation is too slow to avoid extinction, a relative (or quanti- 
tative) constraint is effectively transformed into an absolute 
constraint. They coined the term 'demographic constraint' 
to refer to this situation and calculated 'critical amounts of 
genetic variance' and 'critical heritabilities' that are neces- 
sary to prevent extinction under scenarios of sudden and 
gradual environmental change. 

In the multivariate case, an additional source of genetic 
constraints may arise from genetic correlations. Indeed, 
even if every single trait has positive genetic variance, the 
variance for certain trait combinations may be zero (Dick- 
erson 1955). In this case, the G-matrix is singular (Lande 
1979), that is, at least one of its eigenvectors has a zero 
eigenvalue. If the selection gradient is parallel to such an 
eigenvector, it will produce no effect. Regardless of the 
direction of selection, evolution will be possible only in a 
lower-dimensional subspace of the original phenotype 
space (e.g., along a line in two dimensions or a plane in 
three dimensions). A singular G-matrix might be an 
extreme case (and is difficult to infer statistically). How- 
ever, relative constraints arise in the same way, whenever 
an eigenvalue is positive but small. Using their concept of 
demographic constraints, Gomulkiewicz and Houle (2009) 
calculated critical values for the smallest eigenvalue of G in 
the worst-case scenario that selection acts exactly in the 
direction of the corresponding eigenvector. 

What is the overall role of genetic correlations in con- 
straining the rate of adaptation? - Walsh and Blows (2009) 
argued that strong multivariate constraints (weak variation 
in the direction of selection) might, indeed, be common 
and could explain the frequent observation of slow evolu- 
tionary change despite strong selection on (individually) 
variable traits. To quantify the distribution of genetic varia- 
tion, Kirkpatrick (2009) defined a measure of 'effective 
dimensionality' 

«- = Er' (5) 

where the A; denote the eigenvalues of the G-matrix 
ordered from the largest (Aj) to the smallest (/„). If genetic 
variation is uniformly distributed among the eigenvectors, 
n,i takes it maximal value of n, whereas it is minimal (equal 
to 1) when genetic variation is only present along a single 
axis. A review of empirical estimates of n d suggests that it is 



often (much) smaller than the number of traits considered 
(Kirkpatrick 2009). Thus, genetic variation seems to be 
concentrated around a few dimensions, meaning that the 
ability of populations to respond to arbitrary selection 
pressures may be severely reduced. 

However, an alternative approach by Agrawal and 
Stinchcombe (2009) yields more nuanced results. These 
authors proposed to compare the increase in mean fitness 
in response to a given selection gradient for the full G- 
matrix with the expected response when assuming a (hypo- 
thetical) modified G-matrix in which all off-diagonal 
entries (i.e., all covariances) have been set to zero. Using 
data from empirical estimates of G- (or P)-matrices and 
selection gradients, they found that removing genetic corre- 
lations sometimes increases and sometimes decreases the 
rate of adaptation and that often, the effect is minor. In this 
context, it is worth pointing out that genetic correlations 
do not necessarily decrease the variance in a particular 
direction. For example, adding arbitrary covariances to a 
diagonal G-matrix can only increase genetic variation in 
the direction of the leading eigenvector (Horn and Johnson 
1985, p. 194). 

Theoretical studies have used two approaches to quantify 
genetic constraints (for a review of measures, see Walsh 
and Blows 2009). If the selection gradient is known, adapt- 
ability and constraints should be expressed relative to its 
direction. A sophisticated set of measures was proposed by 
Hansen and Houle (2008), who distinguish 'respondability' 
(the magnitude of overall phenotypic change in response to 
selection in a given direction with unit magnitude), 'evolv- 
ability' (the magnitude of response in the direction of selec- 
tion), 'conditional evolvability' (the magnitude of response 
in selected traits if correlated traits are forced to remain 
constant), and 'autonomy' (the fraction of genetic variation 
in a trait that is independent of potentially constraining 
characters). For cases where the direction of selection is not 
known, several authors have calculated mean rates of adap- 
tation over a distribution of possible selection gradients 
(Hansen and Houle 2008; Kirkpatrick 2009; Chevin 2013). 
When the distribution of selection gradients is uniform, 
genetic correlations have no effect on the mean rate of 
adaptation, because high rates in directions of large varia- 
tion are offset by low rates in directions of small variation 
(Hansen and Houle 2008; Kirkpatrick 2009). When the dis- 
tribution of gradients is not uniform, however, the mean 
rate of adaptation is highest if selection gradients tend to 
coincide with directions of large genetic variation (Chevin 
2013). 

Other determinants of adaptive potential 

We now go on to discuss a broader view of adaptive poten- 
tial and evolvability. Sexual reproduction and genetic 
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recombination have long been hypothesized to facilitate 
adaptation to changing environments (e.g., by bringing 
together alleles on the same genome and reducing the 
effect of clonal interference). For a gradual- change model, 
this was confirmed via simulation by Burger (1999) (see 
also Charlesworth 1993; Waxman and Peck 1999). In par- 
ticular, the increase in genetic variance under directional 
selection (see above) is almost absent in asexual popula- 
tions. 

Several theoretical studies have compared adaptation in 
(sexual) haploid and diploid populations, but the results 
are complex. Haploid populations can be expected to 
evolve faster than diploid populations, because selection is 
more efficient in haploids (Orr and Otto 1994; Otto and 
Gerstein 2008), and this was confirmed experimentally in 
yeast (Gerstein et al. 2011). Nevertheless, haploid popula- 
tions were invaded by diploid strains (Gerstein and Otto 
2011). While in this case, the 'cryptic fitness advantage' was 
attributed to negative frequency-dependent selection, a 
more general advantage to diploidy was proposed by Sellis 
et al. (2011). Using the framework of Fisher's geometric 
model, these authors argued that heterozygote advantage is 
a natural consequence of adaptation in diploids, at least in 
populations that are close to a phenotypic optimum. (The 
reason is that mutations often have smaller phenotypic 
effects in heterozygotes than in homozygotes, such that het- 
erozygotes may have a fitness advantage, while homozyg- 
otes already overshoot the optimum - a probability that 
increases with the number of phenotypic dimensions.) As 
heterozygote advantage favors the maintenance of poly- 
morphism, diploids are expected to have higher levels of 
genetic variation, conferring them an increased adaptive 
potential in case of rapid environmental change. Indeed, 
simulations showed that, in fluctuating environments, dip- 
loid populations maintained higher mean fitness than hap- 
loids, despite a larger standing load (Sellis et al. 2011). 

Again using Fisher's geometric model, Orr (2000) argued 
that evolvability is reduced in complex organisms, because 
mutations are more likely to have negative pleiotropic side 
effects. This 'cost of complexity' can, however, be reduced 
by a modular organization (Welch and Waxman 2003). 
Indeed, several studies have concluded that 'effective 
complexity' is low in many organisms (e.g., Martin and 
Lenormand 2006b; Lourenco et al. 2011). Such low dimen- 
sionality/pleiotropy is predicted to increase the proportion 
of beneficial mutations with large effect, which in turn can 
facilitate adaptation (Gomulkiewicz et al. 2010). Gene-net- 
work models also predict that small network size leads to 
an increased rate of adaptation, faster population recovery 
and higher critical rates of environmental change (Malcom, 
2011a, b). Along similar lines, mutational robustness (i.e., 
the probability for genotypes connected by mutations to 
express the same phenotype) can paradoxically increase the 



adaptive potential of a population by allowing synonymous 
genetic variants to accumulate, thus increasing the muta- 
tional neighborhood of a given phenotype (Gavrilets 1997; 
Fontana and Schuster 1998; Wagner 2008; Draghi et al. 
2010). 

Finally, adaptive potential is likely to be influenced by a 
species' evolutionary history. In particular, species that 
have evolved in variable environments are more likely to 
survive future environmental change than species that have 
long lived under very constant conditions. The idea is not 
only that past fluctuations endow a species with increased 
genetic variation (see above), which has been pretested by 
selection in past environments (Masel 2006; Wagner 2007; 
Hayden et al. 2011), different habitats, or even in another 
species (e.g., introgression) (Rieseberg et al. 2003; Barrett 
and Schluter 2008), but also that the species may have 
evolved increased plasticity and a more flexible genetic 
architecture (Hansen 2006). Indeed, the last two points 
might be related. Several recent models have shown that 
genetic networks that evolved to express plasticity also 
allow for faster genetic adaptation (Espinosa-Soto et al. 
2011; Fierst 2011; Draghi and Whitlock 2012). 

On the other hand, species that have evolved under 
highly stable conditions are expected to be the most sensi- 
tive to environmental change (Overgaard et al. 2011). In 
particular, there is concern that tropical ectoterms might 
be unable to resist increasing temperatures (Janzen 1967; 
Ghalambor et al. 2006; Deutsch et al. 2008; McCain 2009; 
Hoffmann et al. 2012; Urban et al. 2014). Indeed, such 
species are characterized by narrow thermal tolerance 
curves (Amarasekare and Savage 2012) and have narrow 
altitudinal ranges (McCain 2009). If genetic variation in 
the optimal temperature is proportional to the width of 
the thermal tolerance curve (as has been demonstrated for 
Drosophila; Kellermann et al. 2009; Schilthuizen and Kel- 
lermann 2014), they should also have reduced critical rates 
of environmental change (Huey and Kingsolver 1993). 
Quantitative predictions about extinction risk are difficult, 
however, because most studies on thermal tolerances pro- 
vide only relative, not absolute, fitnesses (Deutsch et al. 
2008; Martin and Huey 2008; Bonebrake and Mastrandrea 
2010). 

Adaptation in space 

Real populations are distributed in space, and they can 
react to environmental change by migration in addition to 
genetic evolution and plasticity (Parmesan 2006; Schloss 
et al. 2012). Here, we are not primarily interested in range 
shifts, but instead focus on the effects of gene flow on local 
adaptation in changing environments. 

A natural extension of the gradual-change model dis- 
cussed above considers a shifting environmental gradient, 
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that is, an optimum that changes in both space and time. 
Building on earlier models by Pease et al. (1989), Kirkpa- 
trick and Barton (1997) and Polechova et al. (2009), Dup- 
utie et al. (2012) recently investigated adaptation of 
multiple quantitative traits in response to such a shifting 
gradient. In particular, they addressed how multivariate 
genetic constraints and gene flow alter the adaptive poten- 
tial. While gene flow from maladapted populations can 
potentially constrain local adaptation, it may also promote 
population persistence by enabling the exploitation of lar- 
ger geographic ranges and by spreading favorable alleles 
(Schiffers et al. 2013). Consequently, regardless of the 
number of traits under selection, the critical rate of envi- 
ronmental change is maximized when dispersal is neither 
too weak nor too strong (Alleaume-Benharira et al. 2006; 
Duputie et al. 2012). Population persistence also strongly 
depends on the slope of the spatial gradient. When the gra- 
dient is weak (i.e., the loss of fitness per unit space is small), 
the population remains well adapted over a wide range. 
Conversely, a steep gradient constrains the range. In this 
case, population persistence depends heavily on the geo- 
metric relation of the G-matrix, the shape of the fitness 
landscape and the direction of the spatial gradient. In par- 
ticular, adaptive constraints are minimal whenever the spa- 
tial gradient is colrinear with the direction of weakest 
stabilizing selection and largest genetic variance. Similar to 
the 'flying-kite effect' (Jones et al. 2004), Duputie et al. 
(2012) also found that, when there is indirect selection on 
negatively correlated traits, adaptation in one trait can 
cause another trait to develop a spatial gradient in the 
direction opposite to its optimum. When genetic variances 
are allowed to evolve (as a consequence of selection and 
gene flow, see above), univariate models have shown that 
sufficiently large populations can be perfectly adapted over 
their whole range, albeit at the cost of an increased standing 
load (Barton 2001; Polechova et al. 2009; Bridle et al. 
2010). 

The effect of gene flow on the G-matrix has been studied 
by Guillaume and Whitlock (2007). Using a continent- 
island model, these authors showed that a migration rate of 
about one individual per generation increases the size of G 
by up to 3-fold and may cause its shape and orientation to 
'flip' (albeit only over timescales of several hundred genera- 
tions). These effects are particularly pronounced if other 
factors acting on G, such as the input of mutational vari- 
ance and mutational or selective correlations, are weak. 

The effect of phenotypic plasticity on local adaptation 
and the colonization of new habitats has been studied by 
Chevin and Lande (2011) and Thibert-Plante and Hendry 
(2011). Both studies found that plasticity can facilitate col- 
onization of new habitats, especially if it is expressed after 
migration (i.e., juvenile dispersal). However, no studies to 
date have considered the joint effect of plasticity and 



genetic adaptation in spatially explicit models under envi- 
ronmental change. 

Beyond single species 

Real populations do not evolve in isolation but are embed- 
ded in a network of ecological interactions, and so predic- 
tions of responses to climate change should be made in a 
community context. Several studies have investigated the 
effects of interspecific competition on the rate of adapta- 
tion and the likelihood of evolutionary rescue. Both posi- 
tive and negative effects are possible. The presence of 
competitors can reduce the rate of adaptation in a focal 
species by reducing its population size (and, hence, genetic 
variance or mutational input) and by 'blocking' the access 
to new ecological niches (Johansson 2007; Jones 2008; 
Jones and Gomulkiewicz 2012; Osmond and Mazancourt 
2013). This effect increases the lag load, decreases the criti- 
cal rate of environmental change and can contribute to 
species extinctions. On the other hand, competition may 
also facilitate adaptation if a competitor (or predator) 
'pushes' a focal species in the direction of the new opti- 
mum (Jones 2008; Osmond and Mazancourt 2013). 
Osmond and Mazancourt (2013) argue that both effects 
can be found in recent studies of character displacement 
in Darwin's finches (Grant and Grant 2006). Evolution 
may also be sped up by competitive release if climate 
change causes a competitor to go extinct (Poloczanka 
et al. 2008). 

In the presence of a shifting spatial gradient (see above), 
community evolution depends on the interaction of local 
adaptation and dispersal (de Mazancourt et al. 2008; 
Urban, et al. 2012a, b). De Mazancourt et al. (2008) used 
simulations of a multipatch model to show that species 
often shift their range to new habitats rather than adapting 
to their altered current habitat and that this effect is stron- 
ger in species-rich communities. Urban et al. (2012a) use 
the term 'competitive constraint' to describe the situation 
where a local species is prevented from adapting to a 
changing environment because its habitat is being invaded 
by a competing species already adapted to the new condi- 
tions. This effect is a possible explanation for niche conser- 
vatism (Wiens et al. 2010) during contemporary evolution. 
The opposite effect is also possible; however, local adapta- 
tion of a resident species can prevent the establishment of a 
later-arriving invader (monopolization effect, Urban and 
de Meester 2009; Urban et al. 2012a). And even mal- 
adapted residents can slow range expansions of dispersing 
species into newly available habitats ('boxcar effect': species 
can climb climate gradients only as fast as species further 
up the line; Urban et al. 2012b). In summary, predicting 
community response to environmental change requires 
considering the interactions of two local processes (local 
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community dynamics; local adaptation) and two regional 
processes (immigration from regional species pool and 
immigration from regional genotype pool, Urban et al. 
2012a). 

If we move beyond pairwise interactions, both rapid evo- 
lution and phenotypic plasticity have been shown to con- 
tribute to community stability (e.g., Werner and Peacor 
2003; Yamamichi et al. 2011). Kovach-Orr and Fussmann 
(2013) coined the terms 'evolutionary and plastic rescue' to 
describe situations where this enhanced stability prevents 
species extinctions after an environmental change. Finally, 
evolutionary responses to climate change in complex com- 
munities will not always increase the chances in population 
survival, but may instead lead to 'evolutionary suicide' 
(Ferriere and Legendre 2013). 

Conclusions 

We have reviewed theoretical models of adaptation to 
changing environments, with a focus on evolutionary rates 
of quantitative traits. Unlike models of evolutionary rescue 
by single mutations, the majority of quantitative-genetic 
models consider gradual rather than abrupt environmental 
change. Early models for single traits have introduced the 
concept of a critical rate of environmental change or maxi- 
mal sustainable rate of evolution, beyond which long-term 
persistence is not possible (Lynch and Lande 1993; Burger 
and Lynch 1995). Subsequently, this concept has been 
extended to include multivariate selection (Gomulkiewicz 
and Houle 2009), spatial variation (Duputie et al. 2012) 
and phenotypic plasticity (Chevin et al. 2010). Despite the 
added complexity, it seems unlikely that genetic evolution 
can frequently produce rates of change beyond 0.1 haldanes 
for more than a few generations. Higher observed rates are 
thus likely to be due to phenotypic plasticity, or to be 
accompanied by population decline. Empirical tests of this 
theory are challenging (Gomulkiewicz and Shaw 2013), in 
part due to a strong impact of nonselective stochastic fac- 
tors on observed evolutionary rates, and only one study 
(Gienapp et al. 2013) has attempted to estimate the critical 
rate of change for a natural population (for estimates based 
on physiological models and laboratory data, see Huey and 
Kingsolver 1993 and Willi and Hoffman 2008). We hope 
that, in the future, more such estimates will become avail- 
able from well-studied populations. Another promising 
avenue is experimental evolution under gradually changing 
conditions (Collins 2004; Perron et al. 2008; Lindsey et al. 
2013). 

We have also identified four developing areas that sig- 
nificantly increase the realism of the basic models. These 
include the interactions between phenotypic plasticity 
and genetic evolution, the role of genetic architecture for 
the adaptive potential, adaptation to shifting spatial gra- 



dients and the influence of interspecific interactions on 
rates of adaptation. The former two concern mainly 
internal (organismal) features, whereas the latter two are 
about external (environmental) factors. Further integrat- 
ing these various models promises to significantly 
advance our understanding of species adaptations to 
climate change. 

Acknowledgements 

We thank Reinhard Biirger for helpful discussions and 
Andrew Hendry and three anonymous reviewers for con- 
structive comments on this manuscript. This study was 
partially supported by FWF-Grant P 22581-B17 to MK. 

Literature cited 

Agrawal, A. F. and J. R. Stinchcombe 2009. How much do genetic co- 
variances alter the rate of adaptation?. Proceedings of the Royal Soci- 
ety B: Biological Sciences 276:1 183-1 191. 

Alleaume-Benharira, M., I. R. Pen, and O. Ronce 2006. Geographical 
patterns of adaptation within a species' range: interactions 
between drift and gene flow. Journal of Evolutionary Biology 
19:203-215. 

Alvarez-Castro, J. M., M. Kopp, and J. Hermisson 2009. Effects of epista- 
sis and the evolution of genetic architecture: exact results for a 2-locus 
model. Theoretical Population Biology 75:109—122. 

Amarasekare, P. and V. Savage 2012. A framework for elucidating the 
temperature dependence of fitness. American Naturalist 179:178-191. 

Ancel, L. W. 1999. A quantitative model of the Simpson-Baldwin effect. 
Journal of Theoretical Biology 196:197-209. 

Ancel, L. W. 2000. Undermining the Baldwin expediting effect: does phe- 
notypic plasticity accelerate evolution? Theoretical Population Biology 
58:307-319. 

Arnold, S. J., R. Burger, P. A. Hohenlohe, B. C. Ajie, and A. G. Jones 
2008. Understanding the evolution and stability of the G-matrix. Evo- 
lution 62:2451-2461. 

Baldwin, M. J. 1896. A new factor in evolution. American Naturalist 
30:441-451. 

Barrett, R. and A. Hendry 2012. Evolutionary rescue under environmen- 
tal change? In B. Wong, and U. Candolin, eds. Behavioral Responses 
to a Changing Wold, pp. 216-233. Oxford University Press, Oxford. 

Barrett, R. D. and D. Schluter 2008. Adaptation from standing genetic 
variation. Trends in Ecology and Evolution 23:38—44. 

Barton, N. 2001. Adaptation at the edge of a species' range. In J. Silver- 
town, and J. Antonovics, eds. Integrating Ecology and Evolution in a 
Spatial Context, pp. 365—392. Cambridge University Press, Cam- 
bridge. 

Barton, N. H. and D. Keightly 2002. Understanding quantitative genetic 
variation. Nature Reviews Genetics 3:11—21. 

Barton, N. H. and M. Turelli 1987. Adaptive landscapes, genetic distance 
and the evolution of quantitative characters. Genetical Research Cam- 
bridge 49:157-173. 

Barton, N. H. and M. Turelli 1989. Evolutionary quantitative genetics: 
how little do we know? Annual Review of Genetics 23:337—370. 

Bell, G. 2013. Evolutionary rescue and the limits of adaptation. Philo- 
sophical Transactions of The Royal Society B: Biological Sciences 
368:20120080. 



182 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 169-191 



Kopp and Matuszewski 



Rapid evolution theory 



Berner, D. 2012. How much can the orientation of G's eigenvectors tell 
us about genetic constraints? Evolutionary Ecology 2:1834— 1842. 

Bjorklund, M., E. Ranta, V. Kaitala, L. A. Bach, P. Lundberg, and N. C. 
Stenseth 2009. Quantitative trait evolution and environmental change. 
PLoS ONE 4:e4521. 

Blows, M. W. 2007. Complexity for complexity's sake? Journal of Evolu- 
tionary Biology 20:39^14. 

Bonebrake, T. C. and M. D. Mastrandrea 2010. Tolerance adaptation 
and precipitation changes complicate latitudinal patterns of climate 
change impacts. Proceedings of the National Academy of Sciences 
107:12581-12586. 

Bradshaw, W. E. and C. M. Holzapfel 2006. Evolutionary response to 
rapid climate change. Science 312:1477-1478. 

Bridle, J. R., J. Polechova, M. Kawata, and R. K. Butlin 2010. Why is 
adaptation prevented at ecological margins? new insights from indi- 
vidual-based simulations. Ecology Letters 13:485-494. 

Burger, R. 1999. Evolution of genetic variability and the advantage of sex 
and recombination in a changing environment. Genetics 153:1055— 
1069. 

Burger, R. 2000. The Mathematical Theory of Selection, Recombination, 
and Mutation. Wiley, Chichester. 

Burger, R. and M. Lynch 1995. Evolution and extinction in a chang- 
ing environment: a quantitative genetic analysis. Evolution 49:151- 
163. 

Burger, R. and M. Lynch 1997. Adaptation and extinction in changing 
environments. In R. Bijlsma, and V. Loeschke, eds. Environmental 
Stress, Adaptation and Evolution, pp. 209-239. Birkhauser, Basel, 
Switzerland. 

Burger, R., G. P. Wagner, and F. Stettinger 1989. How much heritable 
variation can be maintained in finite populations by a mutation-selec- 
tion balance? Evolution 43:1748-1766. 

van Buskirk, J. and U. K. Steiner 2009. The fitness costs of developmental 
canalization and plasticity. Journal of Evolutionary Biology 22:852— 
860. 

Charlesworth, B. 1993. Directional selection and the evolution of sex and 
recombination. Genetical Research Cambridge 61:205-224. 

Chesson, P. L. and R. R. Warner 1981. Environmental variability pro- 
motes coexistence in lottery competitive systems. American Naturalist 
117:923-943. 

Chevin, L.-M. 2013. Genetic constraints on adaptation to a changing 

environment. Evolution 67:708—721. 
Chevin, L. M. and F. Hospital 2008. Selective sweep at a quantitative trait 

locus in the presence of background genetic variation. Genetics 

180:1645-1660. 

Chevin, L. M. and R. Lande 2010. When do adaptive plasticity and 
genetic evolution prevent extinction of a density-regulated popula- 
tion? Evolution 64:1143-1150. 

Chevin, L. M. and R. Lande 2011. Adaptation to marginal habitats by 
evolution of increased phenotypic plasticity. Journal of Evolutionary 
Biology 24:1462-1476. 

Chevin, L. M., R. Lande, and G. M. Mace 2010. Adaptation, plasticity, 
and extinction in a changing environment: towards a predictive the- 
ory. PLoS Biology 8:el000357. 

Chevin, L.-M., S. Collins, and F. Lef evre 2013. Phenotypic plasticity and 
evolutionary demographic responses to climate change: taking theory 
out to the field. Functional Ecology 27:967-979. 

Collins, S. 2004. Phenotypic consequences of 1,000 generations of selec- 
tion at elevated C02 in a green alga. Nature 431:566-569. 

Collins, S., J. de Meaux, and C. Acquisti 2007. Adaptive walks toward a 
moving optimum. Genetics 176:1089—1099. 



Crispo, E. 2007. The Baldwin effect and genetic assimilation: revisiting 
two mechanisms of evolutionary change mediated by phenotypic plas- 
ticity. Evolution 61:2469-2479. 

Darimont, C. T., S. M. Carlson, M. T. Kinnison, P. C. Paquet, T. E. 
Reimchen, and C. C. Wilmers 2009. Human predators outpace other 
agents of trait change in the wild. Proceedings of the National 
Academy of Sciences 106:952-954. 

De Jong, G. 1995. Phenotypic plasticity as a product of selection in a var- 
iable environment. American Naturalist 145:493—512. 

Deutsch, C. A., J. J. Tewksbury, R. B. Huey, K. S. Sheldon, C. K. Ghalam- 
bor, D. C. Haak, and P. R. Martin 2008. Impacts of climate warming 
on terrestrial ectotherms across latitude. Proceedings of the National 
Academy of Sciences 105:6668-6672. 

DeWitt, T. J., A. Sih, and D. S. Wilson 1998. Costs and limits of pheno- 
typic plasticity. Trends in Ecology and Evolution 13:77—81. 

Dickerson, G. E. 1955. Genetic slippage in response to selection for mul- 
tiple objectives. Cold Spring Harbor Symposia on Quantitative Biol- 
ogy 20:213-224. 

Draghi, J. A. and M. C. Whitlock 2012. Phenotypic plasticity facili- 
tates mutational variance, genetic variance, and evolvability along 
the major axis of environmental variation. Evolution 66:2891— 
2902. 

Draghi, J. A., T. L. Parsons, G. P. Wagner, and J. B. Plotkin 2010. 
Mutational robustness can facilitate adaptation. Nature 463:353— 
355. 

Duputie, A., F. Massol, I. Chuine, M. Kirkpatrick, and O. Ronce 2012. 
How do genetic correlations affect species range shifts in a changing 
environment? Ecology Letters 15:251-259. 

Espinosa-Soto, C, O. C. Martin, and A. Wagner 2011. Phenotypic plas- 
ticity can facilitate adaptive evolution in gene regulatory circuits. 
BMC Evolutionary Biology 11:5. 

Estes, S. and S. J. Arnold 2007. Resolving the paradox of stasis: models 
with stabilizing selection explain evolutionary divergence on all time- 
scales. American Naturalist 169:227-244. 

Ferriere, R. and S. Legendre 2013. Eco- evolutionary feedbacks, adaptive 
dynamics and evolutionary rescue theory. Philosophical Transactions 
of The Royal Society B: Biological Sciences 368:20120081. 

Fierst, J. L. 2011. A history of phenotypic plasticity accelerates adaptation 
to a new environment. Journal of Evolutionary Biology 24:1992-2001. 

Fisher, R. 1930. The Genetical Theory of Natural Selection. Clarendon 
Press, Oxford. 

Fontana, W. and P. Schuster 1998. Continuity in evolution: on the nat- 
ure of transitions. Science 280:1451—1455. 

Franks, M.J., J.J. Weber, and S.N. Aitken 2014. Evolutionary and plastic 
responses to climate change in terrestrial plant populations. 
Evolutionary Applications 7:42—55. 

Gavrilets, S. 1997. Evolution and speciation on holey adaptive 
landscapes. Trends in Ecology and Evolution 12:307—312. 

Gavrilets, S. and S. Scheiner 1993. The genetics of phenotypic plasticity. 
V. Evolution of reaction norm shape. Journal of Evolutionary Biology 
6:31-48. 

Gavrilets, S., A. Vose, M. Barluenga, W. Salzburger, and A. Meyer 2007. 

Case studies and mathematical models of ecological speciation. 1. 

cichlids in a crater lake. Molecular Ecology 16:2893-2909. 
Gerstein, A. C. and S. P. Otto 2011. Cryptic fitness advantage: diploids 

invade haploid populations despite lacking any apparent advantage as 

measured by standard fitness assays. PLoS ONE 6:e26599. 
Gerstein, A. C, L. A. Cleathero, M. A. Mandegar, and S. P. Otto 2011. 

Haploids adapt faster than diploids across a range of environments. 

Journal of Evolutionary Biology 24:531—540. 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 169-191 



183 



Rapid evolution theory 

Ghalambor, C. K., R. B. Huey, P. R. Martin, J. J. Tewksbury, and G. 
Wang 2006. Are mountain passes higher in the tropics? Janzen's 
hypothesis revisited. Integrative and Comparative Biology 46:5—17. 

Ghalambor, C. K., J. K. McKay, S. P. Carroll, and D. N. Reznick 2007. 
Adaptive versus non-adaptive phenotypic plasticity and the potential 
for contemporary adaptation in new environments. Functional Ecol- 
ogy 21:394-407. 

Gibbs, H. L. and P. R. Grant 2006. Oscillating selection on Darwin's 
finches. Science 313:224—226. 

Gienapp, P., C. Teplitsky, J. S. Alho, and J. Merila 2008. Climate change 
and evolution: disentangling environmental and genetic responses. 
Molecular Ecology 17:167-178. 

Gienapp, P., M. Lof, T. E. Reed, J. McNamara, S. Verhulst, and M. E. 
Visser 2013. Predicting demographically sustainable rates of adapta- 
tion: can great tit breeding time keep pace with climate change? Philo- 
sophical Transactions of The Royal Society B: Biological Sciences 
368:20120289. 

Gingerich, P. D. 1983. Rates of evolution: effects of time and temporal 
scaling. Science 222:159—161. 

Gingerich, P. D. 2009. Rates of evolution. Annual Review of Ecology, 
Evolution, and Systematics 40:657-675. 

Gomulkiewicz, R. and R. D. Holt 1995. When does Evolution by Natural 
Selection Prevent Extinction? Evolution 49:201-207. 

Gomulkiewicz, R. and D. Houle 2009. Demographic and genetic con- 
straints on evolution. American Naturalist 174:E218— 229. 

Gomulkiewicz, R. and R. G. Shaw 2013. Evolutionary rescue beyond the 
models. Philosophical Transactions of The Royal Society B: Biological 
Sciences 368:20120093. 

Gomulkiewicz, R., R. D. Holt, M. Barfield, and S. L. Nuismer 2010. 
Genetics, adaptation, and invasion in harsh environments. Evolution- 
ary Applications 3:97-108. 

Gonzalez, A., O. Ronce, R. Ferriere, and M. Hochberg 2013. Evolution- 
ary rescue: an emerging focus at the intersection of ecology and evolu- 
tion. Philosophical Transactions of The Royal Society B: Biological 
Sciences 368:20120404. 

Grant, P. R. and R. B. Grant 2006. Evolution of character displacement 
in Darwin's finches. Science 313:224—226. 

Guillaume, F. and M. C. Whitlock 2007. Effects of migration on the 
genetic covariance matrix. Evolution 61:2398—2409. 

Hairston, N. G., S. Ellner, and C. M. Kearns 1996. Overlapping genera- 
tions: the storage effect and the maintenance of biotic diversity. In O. 
E. Rhodes, R. K. Chesser, and M. H. Smith, eds. Population Dynamics 
in Ecological Space and Time, p. 109. University of Chicago Press, 
Chicago. 

Hairston, N. G., S. P. Ellner, M. A. Geber, T. Yoshida, and J. A. Fox 
2005. Rapid evolution and the convergence of ecological and evolu- 
tionary time. Ecology Letters 8: 1 114—1 127. 

Haldane, J. 1957. The cost of natural selection. lournal of Genetics 
55:511-524. 

Hansen, T. F. 2006. The evolution of genetic architecture. Annual Review 

of Ecology, Evolution, and Systematics 37:123—157. 
Hansen, T. F. and D. Houle 2008. Measuring and comparing evolvability 

and constraints in multivariate characters. Journal of Evolutionary 

Biology 21:1201-1219. 
Hayden, E. J., E. Ferrada, A. Wagner, and A. Wagner 2011. Cryptic 

genetic variation promotes rapid evolutionary adaptation in an RNA 

enzyme. Nature 474:92-95. 
Hellmann, J. J. and M. Pineda-Krch 2007. Constraints and reinforcement 

on adaptation under climate change: selection of genetically correlated 

traits. Biological Conservation 137:599—609. 



Kopp and Matuszewski 

Hendry, A. P. and M. T. Kinnison 1999. Perspective: the pace of modern 
life: measuring rates of contemporary microevolution. Evolution 
53:1637-1653. 

Hendry, A. P., T. J. Farrugia, and M. T. Kinnison 2008. Human influ- 
ences on rates of phenotypic change in wild animal populations. 
Molecular Ecology 17:20-29. 

Hereford, J., T. F. Hansen, and D. Houle 2004. Comparing strengths 
of directional selection: how strong is strong? Evolution 58:2133- 
2143. 

Hermisson, J. and P. S. Pennings 2005. Soft sweeps: molecular popula- 
tion genetics of adaptation from standing genetic variation. Genetics 
169:2335-2352. 

Hermisson, J. and G. P. Wagner 2004. The population genetic theory of 
hidden variation and genetic robustness. Genetics 168:2271-2284. 

Hill, W. G. 2010. Understanding and using quantitative genetic varia- 
tion. Philosophical Transactions of The Royal Society B: Biological 
Sciences 365:73-85. 

Hoffmann, A. A. and C. M. Sgro 2011. Climate change and evolutionary 
adaptation. Nature 470:479-485. 

Hoffmann, A. A., S. L. Chown, and S. Clusella-Trullas 2012. Upper ther- 
mal limits in terrestrial ectotherms: how constrained are they? Func- 
tional Ecology 27:934-949. 

Holt, R. and R. Gomulkiewicz 1997. The evolution of species' niches: a 
population dynamic perspective. In R. G. Othmer, F. R. Adler, M. A. 
Lewis, and J. C. Dallon, eds. Case Studies in Mathematical Modelling: 
Ecology, Physiology, and Cell Biology, pp. 25—50. Prentice Hall, 
Englewood Cliffs, NJ. 

Horn, R. and C. Johnson 1985. Matrix Analysis. Cambridge University 
Press, Cambridge. 

Houle, D., J. Mezey, and P. Galpern 2002. Interpretation of the 
Results of Common Principal Components Analyses. Evolution 
56:433-440. 

Huey, R. B. and J. G. Kingsolver 1993. Evolution of resistance to high 
temperature in ectotherms. American Naturalist 142:S21-S46. 

Huey, R. B., M. R. Kearney, A. Krockenberger, J. A. Holtum, M. Jess, 
and S. E. Williams 2012. Predicting organismal vulnerability to cli- 
mate warming: roles of behaviour, physiology and adaptation. Philo- 
sophical Transactions of The Royal Society B: Biological Sciences 
367:1665-1679. 

Hunt, G. 2012. Measuring rates of phenotypic evolution and the insepa- 
rability of tempo and mode. Paleobiology 38:351—373. 

Husby, A., M. E. Visser, and L. E. B. Kruuk 201 1. Speeding up microevo- 
lution: the effects of increasing temperature on selection and genetic 
variance in a wild bird population. PLoS Biology 9:el000585. 

Janzen, D. H. 1967. Why mountain passes are higher in the tropics. 
American Naturalist 101:233-249. 

Johansson, J. 2007. Evolutionary responses to environmental changes: 
how does competition affect adaptation? Evolution 62:421—435. 

Johnson, T. and N. Barton 2005. Theoretical models of selection and 
mutation on quantitative traits. Proceedings of the Royal Society of 
London Series B 360:1411-1425. 

Jones, A. G. 2008. A theoretical quantitative genetic study of negative 
ecological interactions and extinction times in changing environ- 
ments. BMC Evolutionary Biology 8:1 19. 

Jones, E. I. and R. Gomulkiewicz 2012. Biotic interactions, rapid evolu- 
tion, and the establishment of introduced species. American Naturalist 
179:28-36. 

Jones, A. G., S. J. Arnold, and R. Burger 2003. Stability of the G-matrix 
in a population experiencing pleiotropic mutation, stabilizing selec- 
tion, and genetic drift. Evolution 57:1747—1760. 



184 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 169-191 



Kopp and Matuszewski 



Rapid evolution theory 



Jones, A. G., S. J. Arnold, and R. Burger 2004. Evolution and stability of 
the G-matrix on a landscape with a moving optimum. Evolution 
58:1636-1654. 

Jones, A. G., S. J. Arnold, and R. Burger 2007. The mutation matrix and 
the evolution of evolvability. Evolution 61:727—745. 

Jones, A. G., R. Burger, S. J. Arnold, P. A. Hohenlohe, and J. C. Uyeda 
2012. The effects of stochastic and episodic movement of the opti- 
mum on the evolution of the g-matrix and the response of the trait 
mean to selection. Journal of Evolutionary Biology 25:2010-2031. 

Keightley, P. D. and A. Eyre- Walker 2007. Joint inference of the distribu- 
tion of fitness effects of deleterious mutations and population demog- 
raphy based on nucleotide polymorphism frequencies. Genetics 
177:2251-2261. 

Kellermann, V., B. van Heerwaarden, C. M. Sgro, and A. A. Hoffmann 
2009. Fundamental evolutionary limits in ecological traits drive Dro- 
sophila species distributions. Science 325:1244—1246. 

Kimura, M. 1983. The Neutral Theory of Molecular Evolution. Cam- 
bridge University Press, Cambridge. 

Kinnison, M. T. and N. G. Hairston 2007. Eco-evolutionary conservation 
biology: contemporary evolution and the dynamics of persistence. 
Functional Ecology 21:444-454. 

Kinnison, M. T. and A. Hendry 2001. The pace of modern life II: from 
rates of contemporary microevolution to pattern and process. Geneti- 
ca 112-113:145-164. 

Kirkpatrick, M. 2009. Patterns of quantitative genetic variation in multi- 
ple dimensions. Genetics 136:271—284. 

Kirkpatrick, M. and N. H. Barton 1997. Evolution of a species' range. 
American Naturalist 150:1-23. 

Kirkpatrick, M. and S. Peischl 2013. Evolutionary rescue by beneficial 
mutations in environments that change in space and time. Philosophi- 
cal Transactions of The Royal Society B: Biological Sciences 
368:20120082. 

Kopp, M. and W. Gabriel 2006. The dynamic effects of an inducible 
defense in the Nicholson-Bailey model. Theoretical Population Biol- 
ogy 70:43-55. 

Kopp, M. and J. Hermisson 2007. Adaptation of a quantitative trait to a 
moving optimum. Genetics 176:715—718. 

Kopp, M. and J. Hermisson 2009a. The genetic basis of phenotypic adap- 
tation I: Fixation of beneficial mutations in the moving optimum 
model. Genetics 182:233-249. 

Kopp, M. and J. Hermisson 2009b. The genetic basis of phenotypic 
adaptation II: the distribution of adaptive substitutions in the moving 
optimum model. Genetics 183:1453—1476. 

Kovach-Orr, C. and G. F. Fussmann 2013. Evolutionary and plastic res- 
cue in multitrophic model communities. Philosophical Transactions 
of The Royal Society B: Biological Sciences 368:20120084. 

Lande, R. 1976a. The maintenance of genetic variability by mutation in a 
polygenic character with linked loci. Genetics Research 26:221—235. 

Lande, R. 1976b. Natural selection and random genetic drift in pheno- 
typic evolution. Evolution 30:314-334. 

Lande, R. 1979. Quantitative genetic analysis of multivariate evolution 
applied to braimbody size allometry. Evolution 33:402-416. 

Lande, R. 2009. Adaptation to an extraordinary environment by evolu- 
tion of phenotypic plasticity and genetic assimilation. Journal of Evo- 
lutionary Biology 22:1435-1446. 

Lande, R. and S. J. Arnold 1983. The measurement of selection on corre- 
lated characters. Evolution 37:1210—1226. 

Lande, R. and S. Shannon 1996. The role of genetic variation in adapta- 
tion and population persistence changing environment. Evolution 
50:434-437. 



Lindsey, H. A., J. Gallie, S. Taylor, and B. Kerr 2013. Evolutionary rescue 
from extinction is contingent on a lower rate of environmental 
change. Nature 494:463-467. 

Lourenco, J., N. Galtier, and S. Glemin 2011. Complexity, pleiotropy, 
and the fitness effect of mutations. Evolution 65:1559—1571. 

Lynch, M. 1988. The rate of polygenic mutation. Genetical Research 
Cambridge 51:137-143. 

Lynch, M. and R. Lande 1993. Evolution and extinction in response to envi- 
ronmental change. In P. Kareiva, J. G. Kingsolver, and R. B. Juey, eds. 
Biotic Interactions and Global Change, pp. 43—53. Sinauer, Sunderland. 

Lynch, M., W. Gabriel, and A. M. Wood 1991. Adaptive and demo- 
graphic responses of plankton populations to environmental change. 
Limnology and Oceanography 36:1301-1312. 

Malcom, J. W. 2011a. Smaller, scale-free gene networks increase quanti- 
tative trait heritability and result in faster population recoveiy. PLoS 
ONE 6:el4645. 

Malcom, J. W. 2011b. Smaller gene networks permit longer persistence 
in fast-changing environments. PLoS ONE 6:el4747. 

Martin, T. L. and R. B. Huey 2008. Why suboptimal is optimal: Jensen's 
inequality and ectotherm thermal preferences. American Naturalist 
171:E102-E118. 

Martin, G. and T. Lenormand 2006a. The fitness effect of mutations in 
stressful environments: a survey in the light of fitness landscape mod- 
els. Evolution 60:2413-2427. 

Martin, G. and T. Lenormand 2006b. A general multivariate extension of 
Fisher's geometrical model and the distribution of mutation fitness 
effects across species. Evolution 60:893—907. 

Martin, G. and T. Lenormand 2008. The distribution of beneficial and 
fixed mutation fitness effects close to an optimum. Genetics 179:907- 
916. 

Martin, G., R. Aguilee, J. Ramsayer, O. Kaltz, and O. Ronce 2013. 
The probability of evolutionary rescue: towards a quantitative com- 
parison between theory and evolution experiments. Philosophical 
Transactions of The Royal Society B: Biological Sciences 
368:20120088. 

Masel, J. 2006. Cryptic genetic variation is enriched for potential adapta- 
tions. Genetics 172:1985-1991. 

Maynard Smith, J. 1976. What determines the rate of evolution? Ameri- 
can Naturalist 110:331-338. 

de Mazancourt, C, E. Johnson, and T. G. Barraclough 2008. Biodiversity 
inhibits species' evolutionary responses to changing environments. 
Ecology Letters 11:380-388. 

McCain, C. M. 2009. Vertebrate range sizes indicate that mountains may 
be 'higher' in the tropics. Ecology Letters 12:550-560. 

Merila, J. 2012. Evolution in response to climate change: in pursuit of 
the missing evidence. BioEssays 24:811—818. 

Merila, J. and A. Hendry 2014. Climate change, adaptation, and 
phenotypic plasticity: the problem and the evidence. Evolutionary 
Applications 7:1-14. 

Mezey, J. G. and D. Houle 2003. Comparing G matrices: are common 
principal components informative? Genetics 165:411—425. 

Moczek, A. P. 2007. Developmental capacitance, genetic accom- 
modation, and adaptive evolution. Evolution and Development 
9:299-305. 

Orr, H. A. 2000. Adaptation and the cost of complexity. Evolution 
54:13-20. 

Orr, H. A. 2005. The genetic theory of adaptation: a brief history. Nature 
Revue Genetics 6:1 19—127. 

Orr, H. A. and S. P. Otto 1994. Does diploidy increase the rate of adapta- 
tion? Genetics 136:1475-1480. 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 169-191 



185 



Rapid evolution theory 

Orr, H. A. and R. L. Unckless 2008. Extinction and the genetics of adap- 
tation. American Naturalist 172:160-169. 

Osmond, M. M. and C. Mazancourt 2013. How competition affects evo- 
lutionary rescue. Philosophical Transactions of The Royal Society B: 
Biological Sciences 368:20120085. 

Otto, S. P. and A. C. Gerstein 2008. The evolution of haploidy and dip- 
loidy. Current Biology 18:R1 121-1 124. 

Overgaard, J., T. N. Kristensen, K. A. Mitchell, and A. A. Hoffmann 
2011. Thermal tolerance in widespread and tropical drosophila spe- 
cies: Does phenotypic plasticity increase with latitude? American Nat- 
uralist 178:S80-S96. 

Paenke, I., B. Sendhoff, and T. J. Kawecki 2007. Influence of plasticity 
and learning on evolution under directional selection. American Nat- 
uralist 170:47-58. 

Pal, C. and I. Miklos 1999. Epigenetic inheritance, genetic assimilation 
and speciation. Journal of Theoretical Biology 200:19— 37. 

Parmesan, C. 2006. Ecological and evolutionary responses to recent cli- 
mate change. Annual Review of Ecology, Evolution, and Systematics 
37:637-669. 

Pease, C. M., R. Lande, and J. J. Bull 1989. A model of population 
growth, dispersal and evolution in a changing environment. Ecology 
70:1657-1664. 

Pelletier, G., D. Garant, and H. A. P. 2009. Eco-evolutionary dynamics. 
Philosophical Transactions of The Royal Society B: Biological Sciences 
364:1483-1489. 

Perron, G. G., A. Gonzalez, and A. Buckling 2008. The rate of environ- 
mental change drives adaptation to an antibiotic sink. Journal of Evo- 
lutionary Biology 21:1 724-1 73 1 . 

Pfennig, D. W., M. A.Wund, E. C. Snell-Rood, T. Cruickshank, C. D. 
Schlichting, and A. P. Moczek 2010. Phenotypic plasticity's impacts on 
diversification and speciation. Trends in Ecology and Evolution 25:459-467. 

Pigliucci, M. and J. Kaplan 2006. Making Sense of Evolution. The Uni- 
versity of Chicago Press, Chicago. 

Pigliucci, M. and C. J. Murren 2003. Genetic assimilation and a possible 
evolutionary paradox: can microevolution sometimes be so fast as to 
pass us by? Evolution 57:1455-1464. 

Pigliucci, M., C. J. Murren, and C. D. Schlichting 2006. Phenotypic plas- 
ticity and evolution by genetic assimilation. Journal of Experimental 
Biology 209:2362-2367. 

Polechova, J., N. Barton, and G. Marion 2009. Species' range: adaptation 
in space and time. American Naturalist 174:E186-E204. 

Poloczanka, E. S., S. J. Hawkins, A. J. Southward, and M. T. Burrows 
2008. Modeling the response of populations of competing species to 
climate change. Ecology 89:3138-3149. 

Price, T. D., A. Qvarnstrdm, and D. E. Irwin 2003. The role of pheno- 
typic plasticity in driving genetic evolution. Proceedings of the Royal 
Society of London Series B 270:1433-1440. 

Reed, T. E., R. S. Waples, D. E. Schindler, J. J. Hard, and M. T. Kinnison 
2010. Phenotypic plasticity and population viability: the importance 
of environmental predictability. Proceedings of the Royal Society B: 
Biological Sciences 277:3391-3400. 

Rettelbach, A., J. Hermisson, U. Dieckmann, and M. Kopp 2011. Effects 
of genetic architecture on the evolution of assortative mating under 
frequency-dependent disruptive selection. Theoretical Population 
Biology 79:82-96. 

Rieseberg, L. H., O. Raymond, D. M. Rosenthal, Z. Lai, K. Livingstone, 
T. Nakazato, J. L. Durphy et al. 2003. Major Ecological Transitions in 
Wild Sunflowers Facilitated by Hybridization. Science 301:1211—1216. 

Saccheri, I. and I. Hanski 2006. Natural selection and population dynam- 
ics. Trends in Ecology and Evolution 21:341—347. 



Kopp and Matuszewski 

Schiffers, K., E. C. Bourne, S. Lavergne, W. Thuiller, and J. M. J. Travis 
2013. Limited evolutionary rescue of locally adapted populations fac- 
ing climate change. Philosophical Transactions of The Royal Society 
B: Biological Sciences 368:20120083. 

Schilthuizen, M. and V. Kellermann 2014. Contemporary climate change 
and terrestrial invertebrates: evolutionary versus plastic changes. Evo- 
lutionary Applications 7:56-67. 

Schloss, C. A., T. A. Nunez, and J. J. Lawler 2012. Dispersal will limit 
ability of mammals to track climate change in the Western Hemi- 
sphere. Proceedings of the National Academy of Sciences 109:8606- 
8611. 

Schluter, D. 1996. Adaptive radiation along genetic lines of least resis- 
tance. Evolution 50:1766-1774. 

Sellis, D., B. J. Callahan, D. A. Petrov, and P. W. Messer 2011. Heterozy- 
gote advantage as a natural consequence of adaptation in diploids. 
Proceedings of the National Academy of Sciences 108:20666-20671. 

Svardal, H., C. Rueffler, and J. Hermisson 2011. Comparing environ- 
mental and genetic variance as adaptive response to fluctuating selec- 
tion. Evolution 65:2492-2513. 

Thibert-Plante, X. and A. P. Hendry 2011. The consequences of pheno- 
typic plasticity for ecological speciation. Journal of Evolutionary Biol- 
ogy 24:326-342. 

Tollrian, R. and C. D. Harvell (eds.) 1999. The Ecology and Evolution of 
Inducible Defenses. Princeton University Press, Princeton. 

Tonsor, S. J. and S. M. Scheiner 2007. Plastic trait integration across a 
C02 gradient in Arabidopsis thaliana. American Naturalist 169:E1 19— 
E140. 

Tufto, J. 2000. The evolution of plasticity and nonplastic spatial and 
temporal adaptations in the presence of imperfect environmental cues. 
American Naturalist 156:121-130. 

Turelli, M. 1984. Heritable genetic variation via mutation-selection bal- 
ance: Lerch's zeta meets the abdominal bristle. Theoretical Population 
Biology 25:138-193. 

Uecker, H. and J. Hermisson 201 1. On the fixation process of a beneficial 
mutation in a variable environment. Genetics 188:915-930. 

Uller, T. and H. Helantera 2011. When are genes 'leaders' or 'followers' 
in evolution? Trends in Ecology and Evolution 26:435—436. 

Urban, M. C. and L. de Meester 2009. Community monopolization: 
local adaptation enhances priority effects in an evolving metacom- 
munity. Proceedings of the Royal Society of London Series B 
276:4129-4138. 

Urban, M. C, L. de Meester, M. Vellend, R. Stocks, and J. Vanoverbeke 
2012a. A crucial step toward realism: responses to climate change 
from an evolving metacommunity perspective. Evolutionary Applica- 
tions 5:154-167. 

Urban, M. C, J. J. Tewksbury, and K. S. Sheldon 2012b. On a collision 
course: competition and dispersal differences create no-analogue 
communities and cause extinctions during climate change. Proceed- 
ings of the Royal Society B: Biological Sciences 279:2072-2080.: Bio- 
logical Sciences. 

Urban, M. C, J.L. Richardson, and N.A. Freidenfelds 2014. Plasticity and 
genetic adaptation mediate amphibian and reptile responses to cli- 
mate change. Evolutionary Applications 7:88—103. 

Uyeda, J. C, T. F. Hansen, S. J. Arnold, and J. Pienaar 2011. The mil- 
lion-year wait for macroevolutionary bursts. Proceedings of the 
National Academy of Sciences 108:15908-15913. 

Via, S. and R. Lande 1985. Genotype -environment interaction and the 
evolution of phenotypic plasticity. Evolution 39:505—522. 

Wagner, A. 2007. Robustness and Evolvability in Living Systems. Prince- 
ton University Press, Princeton. 



186 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 169-191 



Kopp and Matuszewski 



Rapid evolution theory 



Wagner, A. 2008. Robustness and evolvability: a paradox resolved. Pro- 
ceedings of the Royal Society B: Biological Sciences 275:91-100. 

Walsh, B. and M. W. Blows 2009. Abundant genetic variation + strong 
selection = multivariate genetic constraints: a geometric view of adapta- 
tion. Annual Review of Ecology, Evolution, and Systematics 40:41-59. 

Waxman, D. and J. R. Peck 1999. Sex and adaptation in a changing envi- 
ronment. Genetics 153:1041-1053. 

Welch, J. J. and D. Waxman 2003. Modularity and the cost of complex- 
ity. Evolution 57:1723-1734. 

Wennersten, L. and A. Forsman 2012. Population-level consequences of 
polymorphism, plasticity and randomized phenotype switching: a 
review of predictions. Biological Reviews 87:156—767. 

Werner, E. E. and S. D. Peacor 2003. A review of trait-mediated indirect 
interactions in ecological communities. Ecology 84:1083—1100. 

West-Eberhard, M. 2003. Developmental Plasticity and Evolution. 
Oxford University Press, Oxford. 

West-Eberhard, M. J. 2005. Phenotypic accommodation: adaptive inno- 
vation due to developmental plasticity. Journal of Experimental Zool- 
ogy (Mol. Dev. Evol.) 304B:610-618. 

Westley, P. A. H. 2011. What invasive species reveal about the rate and 
form of contemporary phenotypic change in nature. American Natu- 
ralist 177:496-509. 

Wiens, J. J., D. D. Ackerly, A. P. Allen, B. L. Anacker, L. B. Buckley, H. 
V. Cornell, E. I. Damschen et al. 2010. Niche conservatism as an emr- 
ging principle in ecology and conservation biology. Ecology Letters 
13:1310-1324. 

Willi, Y. and A. A. Hoffman 2008. Demographic factors and genetic vari- 
ation influence population persistence under environmental change. 
Journal of Evolutionary Biology 22:124-133. 

Williams, S. E., L. P. Shoo, J. L. Isaac, A. A. Hoffmann, and G. Langham 
2008. Towards an integrated framework for assessing the vulnerability 
of species to climate change. PLoS Biology 5 6:2621-2626. 



Wund, M. A. 2012. Assessing the impacts of phenotypic plasticity on 
evolution. Integrative and Comparative Biology 52:5-15. 

Yamamichi, M., T. Yoshida, and A. Sasaki 2011. Comparing the effects 
or rapid evolution and phenotypic plasticity on predator-prey dynam- 
ics. American Naturalist 178:287-304. 

Yeaman, S. and M. C. Whitlock 201 1. The genetic architecture of adapta- 
tion under migration-selection balance. Evolution 65:1897—1911. 

Appendix 

Appendix 1: Adaptation and extinction in the one-dimensional 
moving-optimum model 

The following is a simplified version of the model by 
Burger and Lynch (1995), which assumes the Gaussian fit- 
ness function 

( (z-°t) 2 \ , x 

^ = Sesp(-4^j (Al) 

with 

e t = kt (A2) 

here, z is the phenotype of an individual, w Zj( its fitness at 
time t, 9 t is the optimal phenotype, which increases linearly 
at rate k, and to 2 measures the width of the fitness land- 
scape (i.e., selection is strong if co 2 is small). B is the 
expected number of offspring (absolute fitness) of a per- 
fectly adapted individual, and hence, In B is the maximal 
population growth rate. 



h 2 = 0.2 h 2 = 0.4 h 2 = 0.6 




Figure A1 The critical rate of phenotypic evolution, /c crit = k crit I a p (eqn A7) expressed in haldanes, for the one-dimensional moving-optimum 
model (A2) (after Burger and Lynch 1995), as a function of the width of the fitness function V s = co 2 + a\ (with a\ = 1), for various values of the 
reproductive potential 6. The top row shows results for three different values of heritability h 2 = a 2 /(a 2 + al). In the bottom row, a 2 has been set to 
the value predicted by the stochastic house-of-cards (SHC) approximation under pure stabilizing selection for three values of the population size N. 
The SHC approximation is given by a 2 = 2V m A/ e /(1 + a. 2 N e /V s ) (Burger and Lynch 1995), where V m is the mutational variance, a 2 is the variance of 
the effect of new mutations, and N e a 2BN I (26-1) (Burger and Lynch 1995) is the effective population size. The figures are for V m = 0.001 and 
a 2 = 0.05. The thin dotted line gives the heritability h 2 associated with a 2 (SHC). 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 169-191 



187 



Rapid evolution theory 



Kopp and Matuszewski 



If the trait z is normally distributed in the population with 

isolu 



mean z t and variance a 2 , the mean absolute fitness at time f is 



B 



ol + V s 



exp 



9tY 

2(<r? + v s ; 



(A3) 



with V s = m 2 + a 2 e describing the effective width of the 
fitness landscape (which is somewhat 'smeared out' by 
the environmental variance a 2 ). Equation (A3) shows 
that the maximal fitness B is reduced by two components 
of genetic load: The standing load (the square root term) 
due to standing genetic variation, and the lag load due to 
the deviation of the mean phenotype from the optimum. 

For constant a 2 a population with mean phenotype z t 
evolves according to Lande's equation Az t = o 2 /] t , where 
the directional selection gradient at time t is 



ft 



d\n w, 

dz t 



Bt-z, 
a 2 , + V. 



(A4) 



P t measures the proportional change in log mean fitness 
per unit change of the mean phenotype. 

As outlined in the main text, the population will reach a 
state of dynamic equilibrium, where it follows the opti- 
mum with a constant lag, which is given by 



kt — z t = k- 



K + V, 



(A5) 



(Burger and Lynch 1995). 

At the same time, the population dynamics are governed 
by N t+i = N t w t (e.g., Gomulkiewicz and Holt 1995) or a 
density- dependent version thereof (e.g., Burger and Lynch 



N 0 = 500 



^0= 1000 



1995). In any case, population survival requires that, given 
the equilibrium lag 8*, the equilibrium mean fitness 
W* > 1. This condition yields the critical rate of environ- 
mental change 



2 In B 



2lnB 



V, 



(A6) 



(Burger 2000), where the approximation is valid for weak 
selection (V s > 20; see Burger 2000). When scaled by the 
phenotypic standard deviation, eqn (A6) gives the critical 
rate of phenotypic evolution 



^crit/ °~» 



(A7) 



in haldanes. Figure Al illustrates the value of K cri , as a func- 
tion of heritability, the width of the fitness landscape and 
the reproductive potential of the population. The rule of 
thumb /v crit < 0.1 (Burger and Lynch 1995) is based on V s 
between 5 and 100, In (5) < 1, and h 2 < 0.5. Burger and 
Lynch (1995) note that genetic drift and fluctuating selec- 
tion might decrease rC crit even further. 

Appendix 2: Adaptation and extinction in the multidimensional 
moving-optimum model 

The multivariate version of model (A2) for n selected traits is 
w M = B exp (-^ (z - e t )'a>- 1 (z - 0 f )) . (A8) 

where the fitness landscape is a bell-shaped 'hill', whose 
orientation and dimensions are determined by the (positive 
semi-definite) covariance matrix co (of size nxn). 



N 0 = 10 000 



N 0 = 100 000 



o 

1 05 
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Figure A2 Critical rates of environmental change k CIlt (top row) and the corresponding rates of phenotypic evolution K- crit (bottom row), under the 
premise that the population maintains a minimal size of A/ crit individuals over t crit generations. The case f crit = oo, W crit = W 0 (where A/ 0 is the initial pop- 
ulation size, which equals the carrying capacity; dashed line) corresponds to the case investigated by Burger and Lynch (1995). The case f crit = 50, 
Wcrit = W 0 (gray line) is given by eqn A14. In the bottom row, solid and dashed gray lines are identical, because )c crit (t crit , A/ 0 ) = K C rit (oo, N 0 ). Parame- 
ters are as in the bottom row of Fig. A1 with 6 = 2. 
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If the phenotype distribution is multivariate normal, the 
mean fitness is 

w t = B- S J det(((o + P)- 1 t U )exp ~ 6,)'{m + P) _1 (z f - 0,)\ 

(A9) 

which, as in the univariate case, is reduced by a standing 
load (the square root term) and a lag load (the exponen- 
tial). 

The mean phenotype evolves according to the multivari- 
ate Lande equation Az t = G/? t , where the multivariate 
selection gradient is 



fi t = {m + V)- 1 (6 t -z t ) 



(A10) 



In the gradual-change scenario, 6, = Vt, the rate and the 
direction of environmental change is described by a speed 
vector k=(fc!, . . . which contains the rates of change in 
the optimum for each trait. As in the univariate case, the 
population will eventually follow the optimum with a con- 
stant lag (assuming the G-matrix is constant): 

5* = (<o + P)G _1 k 

(Jones et al. 2004; Gomulkiewicz and Houle 2009; Chevin 
2013; Jones et al. 2012). Again, the population can persist 
if w* > 1. If the vector k is decomposed into its length and 
its direction, k = ||k||c with c = k/||k||, then this condition 
is satisfied if 



/ 2ln(B)+ln(det(( m + P)- 1 t »)) 
c^G-'(« + P)G-'c 

(Gomulkiewicz and Houle 2009). 

Appendix 3: Fixation of a major mutation 

Assume constant selection and a major mutation increasing 
fitness from 1— s to 1. Let the phenotypic effect of this 
mutation be Sz. For simplicity, we look at the haploid case 
and neglect genetic variation at other loci as well as envi- 
ronmental variance. Then, the phenotypic variance of the 
population at time t is ai t = p t (l — pt)dz 2 , where p t is the 
frequency of the beneficial allele. The per generation change 
in mean phenotype is 



Az, = 8zp,(l -p t ) 



1-S+Spt 



which in haldanes is 
Az, 



P , t = vpt{1 - pt) r^p- t 



(A12) 



(A13) 



independent of 8z. The maximal rate of change, which is 
achieved when p = 1/2, is s/(2— s), which may be large if 
selection is strong (e.g., for s = 1/2, K max = 1/2). 



1.00 
0.75 
0.50 
0.25 
0.00 
-0.25 




An= 1000 



N n = 10 000 



n= 100 
n=N 



» ■■Jiiii LyJijiiliLk lliitiiiif 

HtVtillitai^adi mdlilhilLlrihilfrilii Irtish 



_i_ 



100 



200 300 
f 



400 



500 



1.00 
0.75 
0.50 
0.25 
0.00 
-0.25 




| i 'IP T | 1 



I 



_i_ 



100 



200 300 
t 



400 



500 



Figure A3 Observed generation-to-generation rates of phenotypic change k in haldanes for the entire population (n = N, black line) or based on a 
sample of n = 100 individuals (gray line), for two simulation runs with carrying capacities (and initial population sizes) W 0 = 1000 and W 0 = 10000, 
respectively, and parameters as in Fig. A2. The inset shows the trajectories of the mean phenotype z and the phenotypic optimum z opt . The spike in k 
around generation 40, which partially closes the large initial phenotypic lag (insert), is due to an increase in genetic variance (see main text). Fluctua- 
tions in the black line reflect genetic drift and environmental variance, whereas those in the gray line are largely due to sampling effects. In addition, 
rates measured in haldanes vary due to fluctuations in the phenotypic variance oi (for potential problems of scale, see Hereford et al. 2004; Hansen 
and Houle 2008). 
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Appendix 4: Maximal sustainable rates of phenotypic evolution 
over 'time frames of conservation interest' 

Barrett and Hendry (2012) have argued that, over time 
frames of conservation interest, maximal sustainable rates of 
phenotypic evolution could well exceed the 0.1 haldanes that 
have been proposed by Burger and Lynch (1995) for the 
long-term equilibrium of the moving-optimum model. 
Their point was that environments will not keep changing 
forever and that conservation biology is rather concerned 
with population survival over modest periods of time (e.g., 
50 generations). Here, we attempt to evaluate this claim by 
calculating critical rates of environmental and phenotypic 
change, fc cri t(fcrit> N crit ) and fc cri t(fcrit> Ncrit) > such that the 
population consists of N crit individuals after f crit generations. 
Note that, as we are no longer considering a dynamic equi- 
librium, the two rates fc crit and K crit are no longer equivalent. 
Our analysis is based on a 'quasi -deterministic' approxima- 
tion developed by Burger and Lynch (1995) for studying the 
mean time to extinction. This analysis neglects evolution of 
genetic variance in response to selection, as well as several 
sources of stochasticity (see Appendix 5). Its results can 
therefore only be a first approximation, which, however, 
help to elucidate several principals. 

Consider, first, the case N crit = N 0 , that is, we require 
that the population size does not decline from its initial 
value N 0 over f crit generations. The corresponding critical 
rate of environmental change can be calculated by rearrang- 
ing equation (12a) in Burger and Lynch (1995), which gives 



fccrit ( fcrit ,N 0 ) = fccri, (00 , N 0 ) 



(A14) 

where k clit (<x>, N 0 ) is the critical rate given in eqn (A6) for 
infinite times. While A: crit (f crit , N 0 ) can substantially exceed 







^1 — exp 


■> l crit 







/c cr it(°o, No), the corresponding critical rates of phenotypic 
change are identical, that is, /c crit (f crit , N 0 ) = K C rit(°°> 
N 0 ) = fc C rit( 00 ) N 0 )/(T p . The reason is that both K cri ,(f cri „ N 0 ) 
and fc crit (co, N 0 ) are achieved when the mean absolute fitness 
w = 1 (see eqn A3). This illustrates that eqn (A6) does not 
depend on the assumption of an indefinitely moving opti- 
mum and that the corresponding K crit simply gives the maxi- 
mal rate at which the population can evolve without 
decreasing in size. In other words, rates of phenotypic evolu- 
tion can only exceed K cl . it (oo, N 0 ) if population size declines. 

To study this case, we now allow moderate population 
decline while still requiring the population size to remain 
above a critical threshold N crit < N 0 over f crit generations. 
No analytical solution exists (Burger and Lynch 1995), but 
the critical rates can be estimated numerically by iterating 
eqn (4). As a reduction in population size also entails a 
reduced genetic variance, we followed Burger and Lynch 
(1995) by assuming that er? at any time is given by the sto- 
chastic house-of-cards approximation for the current N. 

Figure A2 shows maximal rates of environmental and 
phenotypic change under the constraint that the popula- 
tion is to maintain a minimal size of N crit = 50 individuals 
over t crit = 50 generations and compares the results to 
those from the Burger and Lynch (1995) framework (no 
reduction in population size over infinite times) and those 
from eqn (A14) (no reduction in population size over 50 
generations). The critical rate of environmental change, 
fc crit , increases substantially when short-term reductions in 
population size are allowed (top row of Fig. A5). In con- 
trast, the differences in the critical rates of phenotypic 
change, K clit are much less pronounced (bottom row), in 
particular for small populations (which might be of highest 
interest for conservation). Even in large populations, the 
relative increase in K crlt rarely exceeds 30% (unless selection 
is extremely strong). Given the large uncertainty in our 
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Figure A4 Distribution of observed generation-to-generation rates of phenotypic change k in haldanes, over 100 simulation runs similar to those in 
Fig. A3. In (A), rates based on the entire population (n = N) or on samples of size n = 100 are shown for various initial population sizes (and carrying 
capacities) W 0 . In (B), W 0 = 10000 was kept constant and only sample size n was varied. Other parameters are as in Fig. A1 . 
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estimates of evolutionary rates (see below), this increase 
appears minor, and we conclude that considering adapta- 
tion over 'time frames of conservation interest' does not 
substantially alter the rule-of-thumb that critical rates are 
typically around 0.1 haldanes. 

Appendix 5: Stochastic fluctuations in evolutionary rates 

The analysis presented in Appendix 4 was based on a deter- 
ministic approximation, which neglects various sources of 
stochasticity (see main text). To illustrate this stochasticity, 
we conducted individual-based simulations as described in 
Burger and Lynch (1995). Two exemplary runs are shown 
in Fig. A3. While the population mean phenotype follows 
the moving optimum, generation-to-generation rates of 
phenotypic change (fc) in haldanes fluctuate as a conse- 
quence of nonselective factors such as genetic drift, envi- 
ronmental variance and fluctuations in the phenotypic 
variance a 1 . Observed fluctuations in k are further ampli- 
fied if only a part of the population is sampled (gray lines 
in Fig. A3), and their range can largely surpass the 0.1 hald- 
anes predicted by Burger and Lynch (1995), see Fig. A4. 
Similarly, in small populations, drift and environmental 
variance alone can induce rates of changes of up to 0.15 
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Figure A5 Expected absolute rate of phenotypic change k between 
generations due to genetic drift and environmental variance, which 
contribute /(nW e ) ar| d 2^/(1 - h 2 )l(%N), respectively (see main 
text), as a function of heritability h 2 and population size N = W e (inset). 
Environmental variance a 2 , refers to the phenotypic variance caused by 
developmental instability and micro-environmental fluctuations. 

haldanes (Fig. A5). Overall, these results cast serious doubt 
on our ability to predict the fate of populations based on 
short-term measures of micro -evolutionary change. 
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